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ABSTRACT 

High-power  free  electron  lasers  (FELs),  capable  of  deployment  aboard  naval  combatants, 
would  place  a  unique  and  significant  demand  upon  the  ship's  electrical  distribution  system.  A 
shipboard  FEL  must  be  power  efficient,  relatively  compact,  and  present  a  minimum  radiation 
hazard  to  nearby  personnel.  The  feasibility  of  deploying  an  FEL  aboard  a  ship  is  analyzed 
from  a  power  system  perspective. 

To  produce  1  MW  of  laser  power,  it  is  determined  that  6.6  MW  of  high-voltage  dc 
power  is  required  to  drive  the  FEL  when  superconductor  accelerator  technology  is  employed 
and  9  MW  is  required  when  conventional  room  temperature  accelerator  technology  is  used. 
The  required  prime  power  electrical  distribution  is  easily  compatible  with  the  gas  turbine 
engineering  plants  of  modern  surface  combatants.  This  distribution  will  add  22  tons  to  the 
ship's  displacement  and  require  22  m3  of  the  ship's  volume  to  implement  Simulation  results 
show  that  the  FEL  would  require  an  undulator  with  only  16  periods  to  produce  1  MW  for  the 
electron  beam  parameters  developed  during  the  power  analysis.  This  FEL  exhibits  a  large 
tolerance  to  electron  beam  quality.  From  a  power  analysis  viewpoint,  FELs  may  become  a 
competitive  technology  for  a  prospective  naval  laser  weapon. 
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I.  INTRODUCTION 

The  free  electron  laser  (FEL)  was  first  proposed  by  J.M.J.  Madey  in  1971  [1]  and  lased 
in  1976  at  Stanford  University.  The  FEL  technology  was  so  promising  that  extensive  funding 
and  resources  were  devoted  to  FEL  research  as  part  of  the  Strategic  Defense  Initiative  (SDI). 
Consequently,  rapid  developments  in  FEL  theory  and  design  were  made  throughout  the  1980s. 
Currently,  the  technology  has  matured  to  where  high-average  power  FELs  are  feasible. 
Stanford  University,  Boeing,  Rocketdyne,  and  the  Continuous  Electron  Beam  Accelerator 
Facility  (CEBAF)  are  among  those  pursuing  kilowatt  level  FELs. 

The  simplified  schematic  of  a  generic  FEL  is  depicted  in  Figure  1. 
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Figure  1  -  Simplified  FEL  diagram 

A  relativistic  electron  beam  is  deflected  periodically  by  an  oscillating  magnetic  field  produced 
by  a  series  of  magnets  called  the  undulator.  The  electron  beam  consists  of  a  series  of 
picosecond  electron  pulses,  called  micropulses,  separated  in  time  with  a  duty  cycle  about  10-4. 
Often  the  micropulses  are  produced  in  trains  several  microseconds  long  at  repetition  rates  of  a 
few  hertz.    One  full  cycle  of  this  repetition  is  called  a  macropulse.   However,  continuous  wave 
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(CW)  operation  of  electron  beams,  where  a  continuous  stream  of  micropulses  are  injected  into 
the  undulator,  can  also  be  used.  Light  is  radiated  from  the  accelerated  electrons,  collected  in  an 
optical  resonator  cavity,  and  then  amplified  by  subsequent  electron  pulses. 

Light  recirculates  within  the  optical  resonator  in  pulses.  The  resonator  cavity  dimensions 
are  aligned  so  that  each  optical  pulse  co-propagates  with  an  electron  micropulse  as  it  travels 
down  the  length  of  the  undulator.  The  micropulse  is  several  hundred  optical  wavelengths  long 
so  that  within  a  length  of  the  micropulse  corresponding  to  an  optical  wavelength,  the  electron 
longitudinal  spatial  distribution  is  essentially  homogeneous. 

FELs  are  an  attractive  light  source  for  several  reasons.  Perhaps  the  most  positive 
attribute  of  FELs  is,  unlike  chemical  or  solid-state  lasers,  that  the  optical  wavelength  can  be 
tuned  over  a  wide  range  of  wavelengths.  The  wavelength  is  a  function  of  electron  beam 
energy.  With  slight  variations  in  the  beam  energy,  the  FEL  wavelength  can  be  tuned  over  a 
range  of  about  an  order  of  magnitude  about  the  center  wavelength.  A  more  complete 
description  of  FEL  theory  is  presented  in  Chapter  II. 

The  FEL  also  has  disadvantages.  The  average  power  of  FELs  has  been  limited  to  only  a 
few  watts  while  the  capital  cost  of  the  FEL,  including  the  undulator  magnets  and  a 
radiofrequency  (rf)  accelerator  system  used  to  produce  and  transport  the  electron  beam,  can 
cost  tens  of  millions  of  dollars.  Electron  accelerators  have  been  unable  to  produce  the  large 
average  currents  required  to  generate  high-average  power.  Technology  will  now  support  a 
kilowatt  scale  FEL  and  a  100  kW  FEL  is  expected  within  5-10  years  [2]. 

Free  electron  lasers  are  of  interest  to  the  Department  of  Defense  due  to  their  potential  as 
high-power  weapons  and  lower  power  countermeasure  devices.  A  shipboard  high-power  FEL 
could  provide  essentially  a  zero  time-of-flight  defense  against  enemy  missiles  and  aircraft.  The 
need  for  such  a  weapon  has  been  previously  established  [3] [4].  Technological  developments 
required  to  build  an  FEL  weapon,  however,  are  probably  more  than  ten  years  away. 

A  shipboard  FEL  would  require  an  electrical  prime  power  source.  The  amount  and  type 
of  power  required  by  an  FEL  would  place  a  specific  demand  on  a  ship's  electrical  power 
distribution  system.  Several  possible  choices  of  FEL  system  architectures  are  potential 
candidates  for  a  high-power  FEL.  The  maritime  combat  environment,  however,  places  unique 
design  constraints  on  the  system.  In  Chapter  III,  the  power  requirements  of  the  FEL 
architecture  most  suited  to  this  environment  are  described.  The  methods  used  to  determine 
these  requirements  are  also  valid  for  other  architectures  with  appropriate  changes  in  design 
parameters.  All  of  the  technologies  considered  are  in  use  on  existing  FELs  or  have  otherwise 
been  demonstrated.  No  attempt  to  predict  future  technologies  has  been  made.    Chapter  IV 


outlines  the  prime  power  considerations  in  producing  and  distributing  this  power  for  the 
shipboard  engineering  plants  that  will  power  naval  combatants  into  the  21sl  century. 

The  power  analysis  of  Chapter  III  makes  several  assumptions  that  would  affect  the  design 
and  operation  of  the  undulator  and  optical  resonator.  In  Chapter  V,  simulations  are  presented 
to  demonstrate  that  an  FEL  can  operate  at  design  levels  under  these  assumptions.  Current 
accelerator  and  optics  technology,  however,  would  not  likely  be  able  to  support  all  of  these 
parameters  but  it  is  important  to  validate  the  feasibility  of  pursuing  such  a  design  or  a 
derivative. 

This  thesis  incorporates  many  of  the  aspects  of  FEL  design.  An  FEL  is  not  comprised  of 
a  set  of  unique  characteristics.  There  are  many  design  choices  which  must  be  made  in  order  to 
define  a  given  FEL  architecture.  The  feasibility  of  a  high-power  shipboard  FEL  weapon  is 
examined  from  a  power  analysis  approach.  A  model  for  an  FEL  weapon  is  defined  based  on 
the  volume,  power  distribution,  and  environmental  constraints  of  the  maritime  combat 
environment.  The  power  requirements  for  this  model  are  determined  and  the  characteristics  of 
the  required  power  distribution  system  are  outlined.  Finally,  the  optical  characteristics  and 
performance  of  this  FEL  model  is  explored. 


II.  FREE  ELECTRON  LASER  THEORY 

The  basic  components  of  an  FEL  are  a  relaiivistic  electron  beam,  a  static  periodic 
magnetic  field,  and  a  co-propagating  optical  field.  The  FEL  operates  through  the  quantum- 
mechanical  processes  of  spontaneous  emission,  stimulated  emission  and  stimulated  absorption. 
The  quantum-mechanical  descriptions  of  these  processes  for  bound  electrons,  as  in  an  atomic 
laser,  are  well-known  and  easily  formulated.  It  is  not  necessary,  however,  to  use  the 
quantum-mechanical  descriptions  for  the  FEL.  The  FEL  can  be  completely  described  by 
classical  electromagnetic  theory  relations  since  there  is  a  high  density  of  electrons  within  an 
FEL  electron  beam  and  photons  within  the  optical  wave  [5].  The  following  theory  has  been 
collected  from  a  variety  of  sources.  The  principal  sources  for  each  of  the  following  sections 
are  noted  in  the  subject  headings.  The  following  relationships  are  presented  in  the  cgs  unit 
system. 

A.  RESONANCE  CONDITION  [5][6] 

An  accelerated  electron  radiates  energy.  In  the  FEL,  the  free  electrons  are  accelerated  by 
the  static  periodic  magnetic  field  within  the  undulator.  The  electrons  are  traveling  down  the 
axis  of  the  undulator  at  relativistic  speeds.  Therefore,  the  period  (wavelength)  of  the  magnetic 
fields  as  they  appear  to  the  electron  are  not  the  same  as  the  apparent  undulator  period  as  seen 
by  an  observer  in  the  laboratory  frame.  This  is  often  referred  to  as  the  relativistic  Doppler 
shift.  Since  the  undulator  wavelength  as  seen  by  the  electrons  is  a  function  of  their  energy,  the 
wavelength  of  the  photons  emitted  by  the  electrons  is  also  a  function  of  their  energy.  We  can 
describe  these  effects,  however,  within  the  laboratory  frame  of  reference  due  to  the  relativistic 
invariance  of  our  classical  relationships  [7]. 

The  electron  beam  interactions  with  the  electromagnetic  fields  within  the  laser  cavity  are 
described  by  the  Lorentz  equations  and  the  wave  equation.  The  relativistic  form  of  the  Lorentz 
equations  are: 

—  (ywcff)    =   -ei^  +  ^xB)  (2.1) 

at 

-fCymc)   =   -e$£  (2.2) 

at 

where  y  is  the  Lorentz  factor,  p  =  v7r.  v*  is  the  electron  velocity,  c  is  the  speed  of  light,  m  is 


the  rest  mass  of  an  electron,  e  is  elementary  charge  magnitude  of  an  electron,  E*  is  the  electric 
field  within  the  laser  cavity,  and  B  is  the  total  of  magnetic  fields  within  the  laser  cavity. 

Equation  2.1  represents  conservation  of  momentum  for  the  system  and  is  commonly 
known  as  the  Lorentz  force  equation.  Equation  2.2  describes  energy  transfer  and  would  be 
called  stimulated  emission  or  absorption  in  a  quantum-mechanical  formulation.  For  an  FEL  we 
desire  a  net  transfer  of  energy  from  the  electron  beam  to  the  optical  field  that  is  initially 
produced  from  spontaneous  emission.  This  energy  transfer  is  many  orders  of  magnitude  larger 
than  the  spontaneous  emission  contribution  per  pass  of  the  electron  beam  through  the 
undulator.  We  must  ensure  that  stimulated  emission  occurs  preferentially  to  stimulated 
absorption  in  order  to  make  the  FEL  a  useful  device. 

The  electric  field,  E,  within  the  undulator  is  due  to  the  presence  of  the  optical  field. 
Equation  2.2  shows  that  for  VE  >  0,  energy  will  be  transferred  from  an  electron  to  the  optical 
field.  Furthermore,  special  relativity  states  that  the  speed  of  an  electron  must  be  less  than  the 
speed  of  light  in  free  space  so  that  the  electron  must  travel  slower  than  ("slips"  along)  the  co- 
propagating  optical  field.  The  transverse  component  of  the  electron  velocity  is  periodic  due  to 
the  acceleration  of  the  electrons  by  the  undulator.  While  traveling  through  one  undulator 
wavelength,  \0,  the  electron's  transverse  velocity,  v.,  will  oscillate  one  full  period  as  shown  in 
Figure  2.  In  this  figure,  an  electron's  longitudinal  travel  is  shown  with  it's  average  velocity 
over  one  undulator  wavelength,  Pyc .  After  one  undulator  period,  the  electron  slips  behind  the 
optical  wave  by  one  optical  wavelength,  X.  At  position  z  =  0,  this  electron  is  located  where 
E  =  0  and  v_j_  =  0.  As  the  electron  moves  down  the  undulator,  Vi  increases  and  E  becomes 
positive  so  that  v*-£  >  0  from  0  <  :  <  X0I2.  As  the  electron  continues  down  the  undulator, 
both  v*^  and  E  are  less  than  zero  so  that  v*  E  >  0  from  \0  /2  <  z  <~k0.  The  condition  where 
the  electron  slips  behind  the  optical  wave  by  one  optical  wavelength  as  the  electron  traverses 
an  undulator  wavelength  is  called  the  resonance  condition.  A  single  resonant  electron  transfers 
energy  from  the  beam  to  the  optical  field. 

The  resonance  condition  defines  an  important  feature  of  the  FEL;  specifically,  its 
wavelength  tunability.  To  illustrate  this,  we  begin  with  the  Lorentz  factor  equation.  Since 
P    p  =  p\r  +  p2,  the  Lorentz  factor  (equation  2.3)  can  be  rearranged  to  solve  for 

p\  =  Vi  -a  +r$;)ir  .  (2.4) 

Define  K2  =  y2^2  ■  undulator  parameter.    Typical  values  of  K  and  y  are  K  ~  1  and  y  ~  50- 
500  so  that  1  +  Y^f/Y2  «  1.   With  a  Taylor  expansion,  we  have 
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Figure  2  -  Electron  slippage  of  an  optical  wavelength  as  it  traverses  a  distance  of 
one  undulator  wavelength.  From  Ref.  [4] 
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The  difference  in  velocity  between  the  optical  wave  and  the  electrons  is  r(l-P7),  while  the 
time  it  takes  for  the  electron  to  traverse  an  undulator  period  is  \0l$2c .  Therefore,  the 
resonance  condition  is 

I  +  K2 
A.  =  K  ljL^  •  (2-6) 

2r 

This  mathematical  representation  of  the  resonance  condition  demonstrates  the  wavelength 
tunability  of  the  FEL.    In  section  II. C,  we  will  see  that 

eB„A.„ 

K  =  -~-~  (2.7) 

2rcmcz 

where  B0  is  the  undulator  rms  magnetic  field  strength.  The  undulator  field  strength  and 
wavelength  characteristics  are  fixed  after  undulator  construction.  Therefore,  since  X  «  1/y2,  the 
FEL  optical  wavelength  can  be  rapidly  adjusted  over  a  wide  range  with  small  adjustments  to 
the  electron  energy.    Actual  FEL  experiments  have  produced  optical  radiation  at  wavelengths 


of  10mm  >  X>  300nm[5]. 

B.  SPONTANEOUS  EMISSION  [5][7] 

An  accelerated  electron  will  radiate  electromagnetic  energy.  The  power  radiated  by  a 
relativistic  electron  is  given  by  the  generalized  Larmor  formula  or  Lienard  result.  For  an 
accelerating  force  perpendicular  to  the  direction  of  electron  motion,  the  total  power  radiated  by 
a  single  electron  is 
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This  power  is  distributed  in  space  according  to 

dPe   =  2e¥         gf 

dQ.  nc     (l+y^2)3 


1  - 


4y282cos2<|) 
U+yW 


(2.9) 


where  6  is  the  angular  deviation  from  the  longitudinal  axis  in  the  plane  of  acceleration,  <f>  is  the 
angular  deviation  from  the  longitudinal  axis  perpendicular  to  the  plane  of  acceleration,  and  Q  is 
the  solid  angle  into  which  the  power  is  radiated.  This  power  distribution  is  illustrated  in  Figure 
3.  The  total  power  is  radiated  into  a  narrow  cone  of  width,  A0  a  1/y  and  solid  angle 
AQ  a  1/y2  centered  along  the  direction  of  motion.  This  narrow  beam  causes  the  optical 
radiation  within  the  FEL  laser  cavity  to  have  an  intense  energy  density  complicating  FEL 
resonator  designs  due  to  the  limitations  of  the  optical  components. 

Photons  emitted  from  the  electron  will  have  an  energy,  £  =  he  Ik,  where  A.  is  given  by 
the  resonance  condition.   Therefore, 

2/idLv2 

E= °Ar  (2.10) 

(1+K2) 

where  k0  is  the  wavenumber  associated  with  the  undulator  static  magnetic  field.  The  time  the 
electron  spends  within  the  undulator  interaction  region  is  Ar  =  Lie,  where  L  is  the  total  length 
of  the  undulator.  The  number  of  photons  emitted  per  electron  in  one  pass  of  the  undulator, 
then,  can  be  easily  determined  by  multiplying  the  total  power  emitted  from  an  electron, 
equation  2.8,  by  At  and  then  dividing  by  the  energy  per  photon,  equation  2.10.  In  the  next 
section  we  will  see  that  (32  =  Ar2A:02c2/v2  so  that  the  number  of  photons  spontaneously  emitted 
per  pass,  We ,  is 

We  =  ^  NaK2(l+K2)  ~2NaK2  (2.11) 
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Figure   3    -    Radiation   pattern   from    a   relativistic   electron   (y  =  100)    under    a 
transverse  acceleration,  (a)  <J>  =  0,  (b)  <J>  =  rc/2 


where  a  =  e2\T\c  ~  1/137  is  the  fine  structure  constant  and  N  =  LIX0  is  the  total  number  of 
undulator  periods.  For  an  undulator  design  with  K  ~  1  and  N  -  100,  about  three  photons  will 
be  emitted  per  electron  per  pass  through  the  undulator.  Not  all  of  these  photons,  however  are 
emitted  into  the  coherent  optical  mode.  The  coherent  optical  mode  size  is  controlled  by  the 
Rayleigh  length,  z0,  which  is  the  propagation  distance  required  for  the  optical  mode  cross- 
sectional  area  to  double.  The  number  of  photons  emitted  into  the  solid  angle  &Q.  «  1/Ny2 
defined  by  the  coherent  optical  mode  with  a  Rayleigh  length  comparable  to  the  length  of  the 
undulator  [5]  is 


Wc  =  —  cxK2  (1+K2)  «  aK 

f       3 


(2.12) 


Therefore,  only  about  1%  of  the  electrons  will  spontaneously  emit  a  photon  into  the  coherent 
optical  mode  per  pass. 

The  above  relations  imply  two  significant  factors.  First,  the  total  photon  energy 
emissions  are  many  orders  of  magnitude  less  than  the  total  electron  energy  so  that  momentum 
recoil  imparted  to  the  electron  has  a  negligible  effect  on  the  electron's  trajectory.  Secondly,  the 


process  of  spontaneous  emission  constitutes  the  initial  mechanism  for  FEL  startup  only.  The 
primary  power  transfer  mechanism  from  the  electrons  to  the  optical  fields  after  a  small  power 
density  has  built  up  within  the  laser  resonator  cavity  is  stimulated  emission.  At  saturation,  up 
to  a  few  percent  of  the  electron  beam  energy  may  be  extracted  per  pass  of  a  micropulse 
through  the  undulator  by  stimulated  emission. 

The  electron  micropulses  passing  through  the  undulator  create  optical  pulses  of  similar 
length.  The  expression  for  the  resonance  condition,  equation  2.6,  is  derived  based  on  the 
electron  slipping  one  optical  wavelength  behind  the  photon  every  undulator  wavelength. 
Therefore,  we  might  expect  that  the  minimum  length  of  the  optical  pulse  will  be  NX.  This 
length  is  called  the  "slippage  distance"  since  the  electron  micropulse  slips  a  total  distance  of 
N  k  along  the  optical  pulse  as  the  micropulse  travels  down  the  undulator.  Actual  optical  pulses 
are  often  longer  than  the  slippage  distance,  and  under  certain  conditions,  called  "short  pulse 
effects,"  may  even  be  shorter  than  the  slippage  distance. 

C.  ELECTRON  DYNAMICS  AND  THE  PENDULUM  EQUATION  [5] 

The  energy  transfer  between  the  electron  beam  and  the  optical  wave  occurs  in  the 
interaction  region  of  the  undulator.  Within  this  region,  the  electric  field  is  due  exclusively  to 
the  optical  electromagnetic  wave  while  the  magnetic  fields  are  due  both  to  the  optical  wave  and 
the  static  undulator  magnetic  field.  The  first  FEL  was  developed  with  helically  polarized 
undulator  magnetic  fields  and  has  a  simpler  mathematical  development  for  some  FEL 
phenomena.  Currently  most  FELs  have  linearly-polarized  undulators.  The  fields  within  an 
undulator  linearly  polarized  in  the  vertical  (v)  direction  are  illustrated  by  Figure  4  and  are 
given  by 

gm  =B0  sin(*0z)   v  (2.13) 

£T  =E  cos  (kz  -cor  -Kj>)  x  (2.14) 

Br  =£  sin(kz-<i)t+§)   y  (2.15) 

where  co  =  kc  is  the  angular  frequency  of  the  optical  radiation,  the  subscript  "m "  represents  an 
undulator  parameter  and  the  subscript  "r "  represents  an  optical  wave  parameter. 

In  order  to  describe  the  electron  dynamics  and  trajectories,  we  return  to  the  Lorentz 
equations  as  applied  to  a  linearly  polarized  undulator. 


at  mc  L 
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Figure  4  -  Fundamental  physics  of  the  FEL  interaction.  After  Ref.  [4] 
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dy  _       e 
dt         mc 

t't 

1 

V  ~~  

1    Vi-U-tf 

Equations  2.16  and  2.18  combine  to  obtain 

Pi    =               cos(k0z) 
Y 

eBo  K 

2iunc l 

(2.17) 
(2.18) 


(2.19) 

(2.20) 

P,  =  Vl  -[1  +K\\  +cos(2*0r)]/y2  (2.21) 

and  substituting  equation  2.19  into  equation  2.17  we  obtain 

dX  m  V2  -^  co5(itz-cor-Ht»)  .  (2.22) 

dt  yrtic 

Integrating  equation  2.19  and  plugging  in  representative  values  for  existing  FELs  [51,  we  find 
that  the  transverse  deflection  of  the  electron  beam  as  it  passes  through  the  undulator  is  on  the 
order  of  0.2  mm.  Furthermore,  undulator  wavelength  values  are  on  the  order  of  5  cm. 
Therefore,  even  with  N  =  100  periods  for  the  undulator,  the  undulator  is  only  a  few  meters 
long.  The  laser  resonator,  however,  may  be  much  longer  to  allow  the  optical  beam  size  to 
increase  and  reduce  the  power  density  impacting  on  the  mirrors. 

The  electron  trajectories  described  by  the  above  equations  have  oscillatory  components 
within  an  undulator  wavelength.  Over  an  undulator  period,  however,  the  electrons  have  an 
apparent  average  or  constant  trajectory.  It  is  convenient  to  define  the  electron's  position  in 
terms  of  the  electron  phase,  £,  as 

C  ■  (*  +  k0  )z  -  of  (2.23) 

where  z  is  the  position  of  the  electron  given  its  average  motion  neglecting  the  oscillations. 
Since  k  »  k0,  A£  ~  k£z.  A  change  Ar  «  X  corresponds  to  a  A(  ~  In  radians.  From  equation 
2.21,  we  obtain 


=  ^-^f 


and 


r 


z  -I  -4  sin(2k0a)  (2.25) 

a: 
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where  %  =  K2/2(l+K2)  and  k0z  ~  k0ct  since  the  electrons  are  traveling  at  nearly  the  speed  of 
light.    Substituting  equations  2.23  and  2.25  into  equation  2.22  we  obtain 

4l  =     *EK     JCos[(C-Hj)Ksin(2A:0c/)]+cos[(C+(|))-2^c/-^sin(2^c/)]L         (2.26) 
dt        V2  ymc   \ 

Averaging  equation  2.26  over  one  undulator  period  results  in 

— L  -  p cos(C+<i>)  (2.27) 

at  V2  yinc 

Differentiating  equation  2.23  with  respect  to  time  twice  and  equation  2.24  once,  they  combine 
with  the  resonance  condition  and  simplify  to  give 

d±  m  o        HV  ^.^  (2  28) 

dt  y~m 

where  KL (E,)  =  K [J0 (^)-J  j(^)].   Equation  2.28  is  called  the  FEL  pendulum  equation. 

The  pendulum  equation  describes  an  important  physical  phenomena  in  FELs  called 
"electron  bunching".  This  is  easily  seen  in  equation  2.27.  Consider  a  group  of  electrons 
uniformly  distributed  over  a  A£  ~  In  space.  Those  electrons  in  the  vicinity  of  £+<t>  ~  re  will 
experience  a  loss  in  energy  and  slow  down.  Those  electrons  near  £-kJ>  ~  0  will  gain  energy 
and  speed  up.  Therefore,  the  electrons  will  tend  to  bunch  together  near  X,  =  rc/2.  Notice  that 
the  electron  bunching  is  periodic  over  a  A£  =  2n  corresponding  to  AT  ~  k.  The  coupling 
between  the  pendulum  equation  and  the  wave  equation  allows  for  the  net  transfer  of  energy 
from  a  bunched  electron  beam. 

D.  THE  FEL  SELF-CONSISTENT  WAVE  EQUATION  [5][8] 

The  propagation  of  the  optical  wave  is  governed  by  the  wave  equation.  The  form  of  the 
wave  equation  applicable  in  the  undulator  interaction  area  is 


c2dt2 


— >  4tt   — > 

A    =   -— 7,  (2.29) 

c      4 


where  A  is  the  vector  potential  and  7^  is  the  transverse  electron  beam  current.    For  an  electron 
beam  and  optical  wave  in  the  interaction  region  of  a  linearly  polarized  undulator 

X  =  £(z'r)  sin(kz-ut+b)   x  (2.30) 

TXi  =-ec$l  b\x-r,)  x  (2.31) 
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P;  = cos(k0z)  x 


(2.32) 


where  53(.x)  is  the  three-dimensional  delta  function  and  the  subscript  "/"  represents  values  of 
the  i'h  electron  such  that  /  =  0  except  at  the  exact  location  of  an  electron.  The  transverse 
beam  current  is  a  summation  of  the  contributions  of  individual  electrons. 

Substituting  equation  2.30  into  equation  2.29  results  in  a  complex  expression.  Since  we 
are  using  the  FEL  to  produce  coherent,  monochromatic  light,  we  can  apply  the  slowly-varying 
amplitude  and  phase  approximation, 

dE 


dz 


dz 


<   k  <|>  , 


dt 

d^ 
dt 


«   co  £  , 


«*:   co  <J> 


(2.33) 


This  approximation  is  useful  for  coherent  monochromatic  sources  because  the  inverse  Fourier 
transform  of  a  narrowband  signal  corresponds  to  a  slowly-varying  envelope  on  the  optical 
wave.   This  approximation  simplifies  the  second-order  wave  equation  to  the  first-order  equation 


dE 

+1 
c 

dE 

dz 

dt 

cosOb  -co/  +§)-E 


dz     c  dt 


In 

s\n{kz  -co/  -hJ>)  = J  i    x 

c 


(2.34) 


Substituting  in  the  expression  for  the  electron  beam  current  and  averaging  over  a  volume  of  the 
electron  beam  that  is  many  optical  wavelengths  long,  2.34  simplifies  to 

-Mk+k^-a*] -il(k-kB)z-a>i]' 


JL+lA. 

dz     c  dt 


Ee">  =  -2<2enK  £ 


+e 


6J(J-r()  (2.35) 


where  [  ~  ]  represents  the  average  over  many  optical  wavelengths.  Using  equations  2.23  and 
2.25  we  note  that  (k+k0)z-(at={,-E,sin(2k0z)  and  (k-k0)z-(Qt=%-2k0z-%sin(2k0z). 
Computing  the  average  in  equation  2.35,  the  wave  equation  is  now 


JL+iJL 
dz     c  dt 


Eei<Sf  =  -2<2neKL{^)  Y—  63(3c-7t) 
,    Y, 


(2.36) 


An  electron  beam  consists  of  a  large  number  of  electrons  making  the  calculation  of  the 
summation  in  2.36  prohibitive.  However,  we  can  take  a  "few"  electrons  over  an  optical 
wavelength  representative  of  the  entire  electron  beam  within  that  wavelength  and  compute  the 
summation.   Therefore,  the  FEL  wave  equation  is 


d_   ±d_ 
dz     c  dt 


Ee1*  =  -2<2neKL  (©   <  p(z  )e~^  I  y  >,=z  ,  c/ 


(2.37) 
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1  p 

where  p  is  the  electron  density  and  <...>=  —  £(...)■ 

Pi-i 

The  wave  equation  and  the  pendulum  equation  are  the  two  important  equations  describing 

the  electron  beam  and  optical  wave  interactions.    However,  it  is  more  convenient  to  non- 

dimensionalize  these  equations  for  FEL  analysis. 

E.  NON-DIMENSIONAL  ANALYSIS  AND  PHASE-SPACE  [5] 

The  pendulum  equation,  2.28,  and  the  wave  equation,  2.37,  are  the  principal  equations  of 
motion  describing  the  electron  beam  and  optical  wave  dynamics  within  an  FEL.  In  most 
instances,  it  is  necessary  to  analyze  FEL  properties  numerically.  Numerical  analysis  is  simpler 
and  results  more  physically  meaningful  when  these  equations  are  used  in  a  non-dimensional 
form. 

A  convenient  choice  for  non-dimensionalizing  the  FEL  equations  is  to  establish  a 
coordinate  system  by  normalizing  time  to  the  time  it  takes  the  optical  wave  (or  electron  beam) 
to  traverse  the  length  of  the  undulator  and  to  reference  distance  according  to  the  position  of  a 
point  traveling  at  the  speed  of  light.  The  transformations  are  x  =  c\L  t  and  z  =  z  -  c  i ,  where 
x  is  the  non-dimensional  time  and  z  is  the  position  with  respect  to  the  reference  position.  The 
coordinate  z  stays  with  a  point  on  the  field  envelope  traveling  at  speed  c.  For  the  special  case 
where  the  electron  pulse  is  many  optical  wavelengths  long  so  that  p(r )  =  p  and  the  electrons 
have  a  uniform  and  constant  energy,  then  the  non-dimensional  pendulum  equation  and  wave 
equation  become 

f4  =  C   =  ,      ,  ««(W)  (2-38) 

oi~  y~mc~ 

-^-(Ee'*)  =  2J2neKL(fy£-    <  e~^  >  (2.39) 

L  ox  y 

Introducing  the  complex  optical  field  envelope  a  =  \a  I  e'*,  and  the  dimensionless 
current,  j,  where  \a\  =  2^2nLNeEKL(^)/-y1mc2  and  j  =  &pN[enLKL(l;)]2rfmc2,  the  non- 
dimensional  pendulum  equation  and  wave  equation  are 

00 

C     =    \a  I  cos(C-HJ))  (2.40) 

a    =   -  j  <  e~'c-  >  (2.41) 


Typically,  the  optical  mode  waist  (radius),  w0  =  yfkzjn,  within  the  undulator  is  larger  than  the 
radius  of  the  electron  beam.    Therefore,  only  part  of  the  optical  wave  is  amplified  by  the 
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electrons.  This  is  compensated  for  by  including  the  "filling  factor,"  F  =  (rb/w0)2,  in  the 
dimensionless  current,  j  — »  jF . 

The  optical  wave  and  electron  evolution  within  the  electron  beam  are  coupled  according 
to  equations  2.40  and  2.41.  The  electron  phase,  X,,  represents  the  electron  position  within  the 
electron  pulse  and  the  change  in  electron  phase,  AX,  »  kAz,  is  indicative  of  the  change  in 
electron  position.  The  electron  phase  is  periodic  over  a  distance  corresponding  to  one  optical 
wavelength.  Therefore,  if  we  analyze  the  X,  evolution  over  a  length  of  one  optical  wavelength 
(2k  radians  in  phase)  of  the  electron  beam,  then  we  know  the  X,  evolution  over  the  entire 
electron  pulse.   It  is  convenient  to  select  the  X,  reference  so  that  <j>  =  0. 

Another  important  parameter  in  our  analysis  of  electron  evolution  is  the  phase  velocity, 

o  _ 

v  =  X,  =  L[(k  +  k0)$,  -  k].  Differentiating  equation  2.24  and  non-dimensionalizing  we  find 
that  Av  =  4nN  Ay/v.  Therefore,  the  phase  velocity  is  indicative  of  the  electron  energy  and  the 
change  in  phase  velocity  is  indicative  of  the  change  in  electron  energy.  A  phase  velocity  of 
v  =  0  corresponds  to  an  electron  exactly  on  resonance,  whereas  an  electron  with  av=n  will 
advance  about  X/2  ahead  of  an  electron  on  resonance  over  the  length  of  the  undulator. 

The  FEL  phase  space  evolution  depicted  in  Figure  5  shows  an  example  of  the  change  in 
electron  energy  and  position  for  an  initial  uniformly  distributed  monoenergetic  electron  beam. 
If  we  can  form  electron  bunching  around  X,  =  re,  then  a  net  amount  of  energy  will  be 
transferred  from  the  electron  beam  to  the  optical  field.  An  electron  bunch  centered  around 
X,  =  0,  however,  will  cause  a  net  energy  to  be  absorbed  by  the  beam  from  the  field.  Figure  5 
shows  the  phase  space  evolution  of  a  resonant  electron  beam,  and  a  beam  slightly  above 
resonance.  A  resonant  electron  beam  has  no  gain.  This  is  because  for  a  uniformly  distributed 
electron  beam  at  resonance,  every  electron  which  contributes  energy  to  the  optical  field  is 
balanced  by  another  electron  which  absorbs  an  equal  amount  of  energy  from  the  field.  An 
electron  slightly  above  resonance  will  contribute  a  net  energy  to  the  optical  field.  This  will  be 
explained  more  fully  in  the  next  section.  Through  the  process  of  mode  competition  in  the 
FEL,  the  optical  wavelength  will  adjust  itself  so  that  the  electron  beam  will  be  slightly  above 
resonance. 

The  "football-shaped"  curve  in  Figure  5  is  called  the  separatrix.  The  separatrix  curve  is 
given  by 

v,2   =   2\a  I  [1  -  sin(C,  +  <t>)]  (2.42) 

where  the  subscript  "s "  refers  to  the  separatrix  values.  An  electron  orbit  enclosed  within  the 
separatrix  is  a  "closed"  orbit.    This  is  equivalent  to  a  pendulum  swinging  side-to-side  with 
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Figure  5  -  Phase  space  evolution  of  an  electron  pulse,  (a)  resonance,  (b)  maximum 
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insufficient  energy  to  reach  the  vertical  position  and  "swing  over  the  top".  An  electron  orbit 
outside  of  the  separatrix  is  an  "open"  orbit.  Since  these  changes  are  periodic,  an  electron 
which  leaves  the  field  of  A£  =  2re  shown  is  replaced  by  another  electron  from  an  adjacent 
A£  =  2rc  section  of  the  micropulse.  This  is  equivalent  to  a  pendulum  swinging  over  the  top  on 
every  revolution  due  to  a  total  energy  greater  than  that  required  to  overcome  the  potential 
energy  of  the  vertical  position.  In  our  FEL,  however,  energy  is  being  transferred  both  to  and 
from  the  electrons  so  they  are  not  necessarily  confined  exclusively  to  one  region  or  the  other. 
As  the  complex  field  strength,  \a  I,  increases,  the  separatrix  becomes  larger  trapping  a  larger 
fraction  of  the  electrons.  For  a  maximum  transfer  of  energy  from  the  electron  beam  to  the 
optical  field,  the  electrons  must  bunch  at  X,  ~  n  in  the  lower  right-hand  section  of  phase  space 
enclosed  within  the  separatrix. 

F.  WEAK  OPTICAL  FIELDS  EVOLUTION  [5] 

From  the  pendulum  and  wave  equations,  it  is  clear  that  electron  bunching  around  £-k}>  =  re 
provides  for  maximum  gain  of  the  optical  fields.  However,  if  the  electron  beam  is  over- 
bunched  such  that  the  bunch  moves  toward  1,+ty  ~  0  in  phase  space,  then  energy  will  be 
absorbed  back  into  the  electron  beam.  The  weak-fields  assumption,  I  a  I  <  re,  is  defined  such 
that  the  electron  beam  is  not  over-bunched.  Over-bunching  is  indicative  of  strong  optical  field 
FEL  interactions. 

Under  weak  fields,  we  can  expand  the  pendulum  equation  and  wave  equation  to  first 
order  in  la  I  and  combine  to  form 


_  U 


i 


a  =  ±L  f    t'F(t')  IaOc-tOI  e     °    d\'  (2.43) 

2   J" 

where 

oo 

fr(x')=jMf(q)e-i^'dq.  (2.44) 

Equation  2.43  is  referred  to  as  the  integral  equation.  F(x')  is  the  characteristic  function  of  the 
beam  quality  distribution  f(q).  The  initial  phase  velocities  of  the  electrons  are  v,  =  v0  +  q 
and  are  distributed  about  v0  according  to  /(<?).  For  a  perfect  beam,  where  all  of  the  electrons 
are  injected  monoenergetically  with  phase  velocity  v0  into  the  undulator  directiy  on  axis  with 
no  transverse  position  or  velocity  components,  then  F  (x)  =  1 . 

The  wave  equation,  2.41,  tells  us  that  for  low  current  conditions,  j  ^  re,  the  optical  field 
does  not  change  significantly  per  pass.  Assuming  perfect  beam  quality,  the  solution  to  the 
integral  equation  gives  an  optical  field  development  per  pass  to  first  order  in  j , 


a (t)  =  a0  \  1  +  j 


2  -  2cos(v0t)  -  v0xsin(v0T) 


2v; 


(2.45) 


and 


4>(x)  =  > 


2sin(v0x)  -  v0x(l  +  cos(v0x)) 


2v; 


+■ 


(2.46) 


where  o0   is  the  initial  complex  field  strength.    The  gain  of  the  optical  field  through  the 
undulator  in  this  low  gain  region  is  given  by 


G(v0>x)  = 


lfl(T)i2-qQ2  _     y 


2-2cos  (v0  t)-v0  x.v//i  (v0  t) 


(2.47) 


The  weak  fields,  low  current  gain  spectrum  is  shown  in  Figure  6. 


****  Gaixi  and  Phase   Curves   **** 
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Figure  6  -  Weak  fields,  low  current  anti-symmetric  gain  spectrum 

The  anti-symmetric  gain  curve  shows  several  important  FEL  properties.  First,  there  is  no  gain 
at  resonance,  v0  =  0.  Secondly,  the  electron  phase  velocity  must  be  slighdy  above  resonance 
to  achieve  gain.  The  maximum  gain  is  G  =  0.135)  and  occurs  at  an  initial  phase  velocity  v0  = 
2.6. 
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Under  high  current,  high  gain  conditions,  ;  »  Tt,  it  can  be  seen  from  the  wave  equation, 
2.41,  that  the  magnitude  of  the  optical  field  changes  appreciably  over  a  single  pass.  The  self- 
consistent  solution  to  the  integral  equation  for  this  region  provides  much  different  results. 
Maintaining  the  assumption  of  perfect  beam  quality,  an  FEL  operating  on  resonance,  v0  =  0, 
will  have  a  complex  field  development  given  by 

|a(T)|  «  ^  e(//2>«Nw2  (2.48) 

3 


4>(T)  = 


1/3 

-  (2.49) 


giving  gain  along  the  undulator  of 


G(T)Sle072)'N3x  (25Q) 


There  is  now  gain  for  an  electron  beam  directly  on  resonance.  Off  resonance,  the  single-pass 
gain  at  the  end  of  the  undulator  is  given  by 

G(vo)  =  l,0/2)1,3V3(1-(v^/3):)  (2.51) 

where  iB  =  (21  j)1  is  the  bunching  time.  The  gain  spectrum  for  the  high  gain  region  is  shown 
in  Figure  7.  As  in  the  low  gain  case,  the  peak  gain  is  still  slightly  above  resonance.  The  phase 
velocity  that  gives  the  peak  gain,  however,  is  now  a  function  of  current.  As  j  ->  °°,  the  gain 
spectrum  approaches  a  symmetric  spectrum  with  the  peak  gain  occuring  at  resonance. 

The  gain  spectrum  bandwidth  can  be  determined  from  equation  2.51.  Defining  the 
bandwidth  as  the  range  of  phase  velocities  where  the  growth  rate  is  reduced  by  10%  [5],  the 
bandwidth  is  Av0  ~  2y1/3.  Using  this  same  definition,  the  gain  bandwidth  for  the  low  gain  case 
is  about  7t/2.  Therefore,  the  bandwidth  for  the  high  gain  case  is  much  larger  than  for  the  low 
gain  case. 

G.  GAIN  DEGRADATION  DUE  TO  BEAM  QUALITY  [5] 

The  previous  gain  analysis  assumed  perfect  beam  quality.  A  real  electron  beam  has  an 
emittance,  e,  where  the  electrons  enter  the  undulator  at  an  angle,  Gr,  with  respect  to  the 
undulator  axis,  or  at  a  transverse  position  slightly  off  axis,  x0 .  The  electrons  will  also  have  a 
small  energy  spread  between  them.  The  emittance  and  energy  spread  of  the  electron  beam 
translates  into  a  spread  in  phase  velocities  in  phase-space.  The  phase  velocity  spread  of  an 
imperfect  electron  beam  due  to  emittance  is  given  by 
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****  Gain  and  Phase  Curvas  **** 
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Figure  7  -  Weak  fields,  high  current  gain  spectrum 


Av  «  --^  (K\\2  +  Y\2)- 

i+at- 


(2.52) 


and  the  phase  velocity  spread  due  to  energy  spread  is  Av  =  4nN  Ay/y. 

Imperfect  beam  quality  reduces  the  single-pass  FEL  gain  by  degrading  electron  bunching. 
Consider  the  integrand  of  the  integral  equation.  The  complex  field  strength  grows  as  the 
electron  beam  travels  down  the  undulator  causing  the  integrand  to  increase.  Therefore,  the  rate 
at  which  the  optical  field  strength  grows  is  faster.  Equation  2.51  and  a  plot  of  equation  2.47 
shows  that  this  feedback  causes  exponential  growth  in  the  FEL  optical  field  strength.  With 
imperfect  beam  quality,  however,  the  F(x)  decays  over  time  so  that  this  feedback  mechanism  is 
attenuated. 

Electrons  enter  the  undulator  with  a  spread  in  phase  velocities  due  to  a  spread  in 
energies.  Typical  energy  spreads  from  an  accelerator  are  on  the  order  of  Ay/y  ~  0.001.  The 
typical  distribution  of  electron  phase  velocities  is  Gaussian  about  v0 , 


fc(Q)  = 


V2n  or. 


(2.53) 
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where  aG  =  4nN  Ay/y  is  the  standard  deviation.   The  corresponding  characteristic  function  is 

Fc(T)  =  e-°at2/2  (2.54) 

FG  (t)  decays  over  the  length  of  the  undulator  with  a  time  on  the  order  of  l/oc . 

The  change  in  phase  velocity  due  to  beam  emirtance  has  a  quadratic  dependence  on 
electron  position  and  entrance  angle.  Therefore,  f  (q)  is  not  the  same  as  the  distribution  of  x0 
and  Qj. .  Electron  pulses  normally  have  a  Gaussian  distribution  in  space.  A  beam  with  matched 
contributions  from  position  and  angle,  y0.r  =  Kk0x0,  where  Qx  and  x0  are  the  rms  deviations  of 
the  electron  beam  entrance  angle  and  position,  will  have  an  exponential  phase  velocity 
distribution  given  by 

fQ{q)  =  ^-eq,a«  q<0  (2.55) 

where  fo(q)  =  0  for  q    >  0  and   oQ  =  4nNy1Q^Kl-k-K2).    The  corresponding  characteristic 
function  is 

Fq(t)=—- - .  (2.56) 

1  -  io0i 

Fq(x)  decays  over  the  length  of  the  undulator  with  a  time  on  the  order  of  l/o0.    A  a0  or  aG  of 
about  n  will  cause  a  random  phase  spread  of  about  n. 

The  effect  of  beam  quality  on  gain  is  illustrated  in  Figure  8  for  the  case  of  electron 
energy  spread.  At  oG  =  0,  the  gain  spectrum  is  identical  to  Figure  6.  As  oc  increases, 
however,  the  gain  spectrum  broadens,  the  phase  velocity  for  peak  gain  increases,  and  most 
importantly  the  peak  gain  is  reduced.   Electron  beam  emirtance  has  a  similar  effects. 

H.  STRONG  OPTICAL  FIELDS  EVOLUTION  AND  SATURATION  [5] 

When  an  FEL  reaches  the  strong  field  regime,  the  gain  per  pass  is  reduced.  The  FEL 
reaches  saturation  when  the  gain  decreases  to  where  it  equals  the  optical  resonator  loss  rate 
including  the  outcoupling  of  the  mirrors.  Strong  fields  can  be  detected  in  phase  space  by  over- 
bunching  of  the  electrons.   This  occurs  when  \a\  >  it. 

In  strong  fields,  the  separatrix  now  traps  a  significant  amount  of  electrons  in  closed 
orbits.  An  electron  bunch  is  quickly  formed  around  X,  +  <J>  ~  n  where  these  electrons  amplify 
the  optical  wave.  As  the  electrons  lose  energy,  their  phase  slips  towards  the  £  +  <J>  =  0  position 
and  energy  is  absorbed  back  from  the  optical  wave.  This  reduces  the  net  gain  per  pass  of  the 
FEL. 
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Figure  8   -   Degradation  of  the   weak   fields,  low  current  anti-symmetric  gain 
spectrum  due  to  electron  beam  energy  spread 


The  gain  degradation  due  to  strong  fields  is  illustrated  in  Figure  9.  As  the  complex  field 
strength  and  height  of  the  separatrix  increases,  so  does  the  gain  degradation.  Careful 
examination  of  Figure  9  reveals  that  strong-fields  produce  similar  gain  degradation 
characteristics  as  the  energy  spread  effects  of  Figure  8.  This  is  because  the  FEL  interaction 
induces  an  energy  spread  among  the  electrons.  Therefore,  the  gain  is  also  reduced  in  a  manner 
similar  to  a  phase  velocity  spread  of  Av  =  4nN  Av/y.  As  the  total  energy  loss  of  the  electrons 
increases,  a  larger  phase  velocity  spread  is  induced  on  the  electrons.  Eventually  many  electrons 
will  move  near  the  ends  or  even  outside  of  the  gain  spectrum  bandwidth  reducing  the  gain  of 
the  optical  fields.  As  the  optical  fields  approach  saturation,  more  total  energy  is  still  being 
removed  from  each  passing  micropulse,  but  the  fractional  increase  in  optical  energy  (gain)  is 
being  reduced. 


23 


j=7      N=16 

aG=0.4 


Figure  9  -  Gain  degradation  due  to  strong  optical  fields 
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III.  FEL  SYSTEM  POWER  REQUIREMENTS 

The  principles  of  FEL  operation  outlined  in  the  previous  chapter  describe  how  energy  is 
extracted  from  a  relativistic  electron  beam  and  radiated  as  coherent  optical  light.  The 
production  of  a  high  quality  relativistic  electron  beam,  however,  is  a  necessary  engineering 
requirement.  The  basic  functional  components  for  a  high-power  FEL  are  the  electron  injector, 
the  linear  accelerator,  the  undulator  and  optical  resonator,  and  the  electron  beam  dump.  Many 
technologies  and  configurations  for  each  of  these  elements  have  been  developed  and  applied. 
Existing  FELs,  however,  operate  at  only  a  few  watts  of  optical  output  and  do  not  suffer  from 
the  infrastructure  constraints  of  the  naval  shipboard  environment  including  space  availability, 
power  generation  and  distribution,  and  adverse  environmental  conditions.  In  this  chapter,  the 
power  requirements  for  a  high-power  FEL  as  applied  to  a  potential  shipboard  weapons  system 
are  explored.  All  of  the  following  technologies  are  currently  in  use  or  under  development  for 
FEL  applications.  Some  of  these  technologies,  however,  will  need  further  engineering 
development  to  support  high-power  operation. 

A.  SHIPBOARD  ARCHITECTURE 

In  developing  a  shipboard  high-power  FEL,  the  issues  that  are  of  principal  concern  are 
system  power  efficiency,  system  size  and  weight,  and  personnel  radiation  hazards.  The  system 
size  and  weight  are  a  direct  function  of  beam  energy,  linear  accelerator  rf  field  frequency, 
system  power  efficiency,  and  optical  resonator  power  density.  Electron  beam  energy  and  the  rf 
frequency  determine  the  linear  accelerator  size  and  weight  while  the  system  power  efficiency 
and  the  rf  frequency  determine  the  rf  power  source  size  and  weight  characteristics.  Klystrons, 
klystrodes,  and  other  high-power  rf  power  sources  suitable  for  accelerator  applications 
contribute  a  significant  fraction  of  the  system  size  and  weight.  The  optical  resonator  power 
density  is  much  higher  in  an  FEL  than  for  a  chemical  or  solid-state  laser  of  equal  energy 
storage  due  to  the  very  small  beam  size.  Therefore,  FEL  resonator  length  tends  to  be  large  in 
order  to  reduce  the  power  density  at  the  mirrors. 

Of  the  existing  or  proposed,  FEL  architectures,  only  the  energy  recovery  system  depicted 
in  Figure  10  optimizes  all  three  concerns.  In  this  scenario,  electron  pulses  are  injected  into  the 
accelerator  at  some  small  initial  energy  of  approximately  4  MeV.  The  accelerator  increases  the 
beam  energy  by  an  order  of  magnitude  or  more.  Bending  magnets  guide  the  beam  and  direct  it 
into  the  undulator.  The  undulator  will  remove  approximately  2  percent  of  the  beam  energy  and 
convert  it  into  optical  energy.  Additional  bending  magnets  then  guide  the  beam  back  into  the 
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Figure  10  -  Generic  energy  recovery  FEL  architecture 

accelerator  where  it  enters  180  degrees  out  of  phase  with  the  accelerating  fields.  The 
decelerating  electron  beam  then  transfers  it's  energy  back  into  the  rf  fields.  The  energy 
recovered  from  the  decelerated  electron  beam  is  then  used  to  accelerate  another  electron  pulse 
to  repeat  the  cycle.  The  decelerated  electron  pulses  are  then  directed  into  a  beam  dump  where 
their  residual  energy  is  dissipated.  This  residual  energy  is  about  5  MeV  and  is  below  the 
threshold  of  generating  neutron  radiation  when  it  is  dissipated  in  the  beam  dump. 

This  architecture  is  the  most  appropriate  for  several  reasons.  First,  the  energy  efficiency 
is  increased  dramatically  and  the  physical  size  of  the  system  is  reduced  to  that  required  in  a 
shipboard  environment.  Most  importantly,  it  is  the  only  architecture  which  would  reduce  the 
personnel  radiation  hazards  to  manageable  levels  since  there  will  be  no  neutrons  and  much 
lower  energy  gamma  radiation.  Table  1  is  a  fisting  of  the  neutron  production  threshold  for 
several  elements.  These  elements  are  commonly  used  as  structural  materials  in  shipboard 
applications  or  have  significance  with  respect  to  neutron  production  in  accelerator  beam  dumps. 
Only  beryllium  (Be)  has  a  neutron  production  threshold  below  5.9  MeV.   Therefore,  as  long  as 
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Element 

Threshold 

Element 

Threshold 

Element 

Threshold 

(MeV) 

(MeV) 

(MeV) 

Be 

1.7 

Cu 

9.9 

Mo 

7.4 

Os 

5.9 

Ni 

11.4 

Mn 

10.2 

*Li 

7.3 

Fe 

11.2 

Co 

10.5 

*c 

18.7 

Zn 

10.2 

Cr 

12.0 

Table  1  -  neutron  production  thresholds  [9].  All  elements  that  have  isotopes  of 
>  10%  abundance  with  thresholds  below  6  MeV  are  included.  *  denotes 
elements  with  isotopes  of  less  than  10%  abundance  with  thresholds  below 
6  MeV 


beryllium  is  excluded  from  the  beam  dump  structural  materials,  neutron  generation  is  not  an 
issue.  An  architecture  which  does  not  utilize  energy  recovery  will  contain  an  inherent  large 
radiation  problem  that  will  be  much  more  difficult  to  solve  due  to  the  large  radiation  energies 
involved. 

The  feedback  stability  of  this  type  of  system  has  been  numerically  modeled  and  shown  to 
be  stable  for  several  designs  [10].  However,  more  comprehensive  simulations,  experiments, 
and  possibly  control  systems  engineering  will  be  required  to  ensure  stability  at  the  higher 
currents  examined  in  this  thesis. 

Since  we  are  concerned  with  the  FEL's  potential  as  a  naval  laser  weapon,  the  system 
must  produce  high  output  power  at  a  wavelength  suitable  for  propagation  in  a  maritime 
atmospheric  environment.  The  following  power  analysis  will  be  developed  for  an  FEL  with  a  1 
MW  optical  output  power  at  a  wavelength  of  3  |im.  The  FEL  wavelength  can  be  tuned  over  a 
range  from  this  operating  point  so  we  will  still  have  some  flexibility  in  the  actual  operating 
wavelength.  To  establish  the  operational  requirements,  it  is  assumed  that  this  laser  will  only 
need  to  operate  for  three  shots  of  15  second  duration  within  one  minute  and  ten  shots  over  a 
twenty  minute  interval.  This  requirement  has  no  basis  in  the  time  required  to  destroy  a  missile 
with  a  laser.  Rather,  it  was  chosen  because  a  missile  traveling  at  a  speed  of  0.9  Mach  could  be 
engaged  over  a  closing  distance  of  about  5  km. 

To  estimate  the  system  power  requirements,  it  is  necessary  to  determine  the  electron 
beam   energy    and   current.     The   resonance   condition   of  equation   2.6   gives   the  required 
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relationship  between  the  electron  beam  energy,  the  undulator  spatial  wavelength,  and  the 
optical  wavelength.  For  given  undulator  and  optical  wavelengths,  the  electron  beam  Lorentz 
factor  must  be 


/(l  +K2)  k0 


2k 
Current  and  planned  FELs  are  summarized  in  Table  2  [11].    The  undulator  parameter,  K,  is 


typically  on  the  order  of  unity  so  that  y  ~  *Jk0  Ik. 

The  power  in  an  electron  beam  is  the  product  of  the  electron  beam  energy,  Wb  =  ymc2, 
and  the  average  electron  beam  current,  Iavg .  For  the  architecture  of  Figure  10,  only  about  2 
percent  of  the  beam  energy  can  be  extracted  by  the  undulator  to  allow  for  efficient  energy 
recovery  [3]  [12]  [13].  However,  only  a  few  experiments  have  been  performed  which 
demonstrate  energy  recovery  but  none  of  these  have  determined  the  accelerator  energy  recovery 
efficiency,  x\rec ,  dependence  on  the  undulator  energy  extraction  efficiency,  r|0 .  Therefore,  the 
assumption  that  no  more  than  2%  of  the  beam  energy  can  be  extracted  within  an  energy 
recovery  architecture  may  be  conservative.  Regardless  of  T)0,  an  electron  beam  cannot  be 
decelerated  much  below  3  MeV  because  the  electron  velocity  and  rf  phase  velocity  are  no 
longer  essentially  the  same.  The  maximum  energy  recovery  efficiency  is  given  by 
Tlmax  =  1  -  (3  Mev  )/Wh .  However,  energy  extraction  in  the  wiggler  introduces  an  energy 
(velocity)  spread  among  the  electrons  in  the  beam  which  further  reduces  T\rec.  More 
experiments  are  needed  to  determine  the  dependence  of  T\rec  on  r|0 . 

For  a  1  MW  optical  output  with  r|0  =  2%,  50  MW  of  electron  beam  power  is  required. 
In  order  to  limit  the  undulator  and  optical  resonator  size  and  the  required  beam  current,  a 
maximum  beam  energy  is  desired.  From  equation  3.1,  Wb  =  100  MeV  and  Iavg  =  500  mA  are 
chosen  as  the  design  parameters  for  the  production  of  50  MW  of  beam  power.  For  this  beam 
energy,  a  r\rec  =  95%  is  assumed. 

It  is  desirable  to  maximize  r|0  so  that  the  total  electron  beam  power  can  be  reduced. 
This  will  provide  several  advantages.  First,  by  reducing  the  required  lavg ,  a  high-power  FEL  is 
closer  to  realization.  Current  technologies  can  only  expect  to  produce  average  currents  of  10- 
100  mA  [2] [14].  Second,  reducing  Wh  reduces  the  size,  weight,  and  expense  of  the  accelerator 
structures  and  allows  for  the  production  of  longer  wavelength  light  that  may  be  more 
advantageous  in  a  maritime  environment. 
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FEL 

X 

O- 

Ea 

/ 

N 

Ao 

K 

i[im) 

(MeV) 

(A) 

(cm) 

EXISTING  FELs: 

338 

25(,ls 

6 

2 

42 

7.1 

0.71 

UCSB(mm  FEL) 

StanfordCFIRFEL) 

80-200 

15ps 

4 

8 

50 

1 

0.7 

UCSB(FIR  FEL) 

63 

25M.S 

6 

2 

150 

2.0 

0.13 

Tokyo(UT-FEL) 

43 

lOps 

13 

22 

40 

4 

0.7 

HollandCFELIX) 

40 

5ps 

25 

50 

38 

6.5 

1.5 

Osaka(ISIR) 

40 

30ps 

17 

50 

32 

6 

1 

Bruyeres(ELSA) 

20 

30ps 

18 

100 

30 

3.2 

0.8 

Frascati(LISA) 

15 

7ps 

25 

5 

50 

4.4 

1 

Grumman(CIRFEL) 

14 

5ps 

14 

150 

73 

1.36 

0.2 

Beijing(IHEP) 

10 

4ps 

30 

14 

50 

3 

1 

Orsay(CLIO) 

8 

0.3ps 

50 

80 

48 

4 

1 

LANL(AFEL) 

4-6 

lOps 

15 

200 

24 

1 

0.27 

Darmstadt(IR-FEL) 

5 

2ps 

40 

2.7 

80 

3.2 

1 

Stanford(SCAFEL) 

5 

0.7ps 

37 

10 

72 

3.1 

0.83 

VanderbiltCFELI) 

3 

lps 

43 

50 

47 

2.3 

1 

Duke(Markm) 

3 

3ps 

44 

20 

47 

2.3 

1 

BNL(ATF) 

0.5 

6ps 

50 

100 

70 

0.88 

0.35 

BNL(ATF-UV) 

0.25 

6ps 

70 

100 

70 

0.88 

0.35 

LANL(APEX) 

0.37 

lOps 

46 

135 

73 

1.36 

0.58 

Tsukuba(NIJI-IV) 

0.35 

160ps 

300 

5 

2x42 

7.2 

2 

Orsay(Super-ACO) 

0.35 

250ps 

800 

.1 

2X10 

13 

4 

Okazaki(UVSOR) 

0.3 

6ps 

500 

5 

2x8 

11 

2 

Novosibirsk(VEPP) 

0.24 

35ps 

350 

6 

2x33 

10 

1.6 

PROPOSED  FELs: 

231-600 

CW 

1.7 

0.2 

185 

0.8 

.15 

Florida(CREOL) 

NetherlandsCTEUFEL) 

180 

20ps 

6 

350 

50 

2.5 

1 

Rutgers(IRFEL) 

140 

25ps 

38 

1.4 

50 

20 

1 

Moscow  (Lebedev ) 

100 

20ps 

30 

0.25 

35 

3.2 

0.75 

Tokai(SCARLET) 

40 

40ps 

15 

10 

62 

3.3 

1 

Stanford(FIREFLY) 

40 

2ps 

20 

6 

25 

6 

1 

LBL(IRFEL) 

3-50 

33ps 

55 

60 

40 

5 

1 

CEBAF(IRFEL) 

2-24 

2ps 

50 

100 

25 

6 

3.1 

Boeing(APLE) 

10 

60ps 

17 

140 

101 

2.4 

0.2 

Boeing(APLE) 

10 

18ps 

34 

450 

257 

3.9 

1.2 

StanfordCFEL) 

10 

4ps 

24 

25 

52 

2.6 

0.87 

Osaka(ILT) 

10 

2ps 

9 

100 

30 

.66 

0.3 

UCLA(IRFEL) 

10 

2ps 

20 

200 

40 

1.5 

1 

Novosibirsk  (RTM) 

7 

50ps 

51 

100 

4x40 

9 

1-2 

BNL(HGHG) 

3.35 

lOps 

30 

110 

83 

1.8 

1.4 

CEBAF(UVFEL) 

0.15-2 

0.4ps 

200 

200 

48 

3 

1.5 

Osaka(FELI) 

1 

2ps 

170 

100 

50 

6 

1.26 

Rocketdyne 

0.84 

3ps 

90 

500 

160 

2.4 

1.4 

Dortmund(DELTA) 

0.4 

50ps 

500 

90 

17 

25 

3.1 

ENEA-Frascati 

0.2-0.3 

15ps 

2.3 

4 

8 

2.5 

1 

HarimaCHTD 

0.28 

lOOps 

500 

3 

170 

1.8 

4.2 

BNL(DUVFEL) 

0.075 

6ps 

310 

300 

682 

2.2 

1.54 

Duke(Ring) 

0.05 

lOps 

1000 

350 

2x33 

10 

1.7 

SLAC(LCLS) 

0.0004 

O.lps 

7000 

2500 

723 

8.3 

4.4 

Table  2  -  Summary  of  existing  and  proposed  FELs  [11] 
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B.  ACCELERATOR  STRUCTURE 

A  typical  linear  accelerator  is  comprised  of  a  number  of  rf  cavity  cells  like  those  depicted 
in  Figure  1 1 . 


Superconducting 

B-Faclory  Cell 

Shape 


Copper 

Storage  Ring 

Cell  Shape 


(a) 


(b) 


Figure  1 1  -  Typical  accelerator  cells,  (a)  superconducting  (b)  room  temperature. 
FromRef.  [15] 


Linear  accelerator  structures  are  fabricated  from  both  room  temperature  (RT)  and 
superconducting  (SC)  materials.  SC  structures  reduce  operating  power  requirements  due  to  the 
low  ohmic  losses  in  the  accelerating  cavity  walls.  They  also  have  higher  allowed  operating 
fields  and  accelerating  gradients  as  well  as  much  larger  apertures  between  cells  in  the  cavities. 
An  accelerating  electron  forms  an  image  charge  within  the  accelerator  structure.  As  it  passes 
through  an  aperture  connecting  two  cells,  the  image  charge  is  now  much  closer  to  the  electron 
beam  axis.  These  image  charges  have  a  disruptive  effect  on  the  beam  quality.  The  larger 
apertures  mean  that  image  charges  will  have  a  lesser  disruptive  effect  in  SC  structures  than  RT 
structures.  The  higher  gradients  mean  that  SC  structures  can  accelerate  the  electron  beam  to 
the  desired  Wh  over  shorter  accelerator  lengths  than  RT  structures  reducing  the  total  weight 
and  volume  of  the  structure.  RT  structures,  however,  are  forced  to  have  geometries  that 
intensify  the  harmful  beam  quality  effects  in  an  effort  to  minimize  power  dissipation  [16].  RT 
structures  still  have  applications  because  they  are  less  expensive  and  easier  to  fabricate  than  the 
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SC  structures  and  they  don't  require  a  large  and  inefficient  liquid  helium  refrigeration  system. 

The  development,  theory,  and  design  behind  accelerators  is  complicated  and  requires 
extensive  numerical  modeling.  In  performing  a  power  requirements  analysis,  however, 
accelerating  cavities  can  be  characterized  by  a  few  parameters  for  most  purposes  [16]  [17]  [18]. 
The  significant  parameters  are  the  if  frequency,  /  (MHz),  the  cavity  quality  factor,  Q ,  the 
shunt  impedance,  Z  =  Rsh  \Q  (Q/m )  where  Rsh  is  the  shunt  resistance,  the  cavity  active  length, 
La  (m),  and  the  average  acceleration  gradient,  Ea  (MV/m).  The  power  dissipated  per  unit 
length  is 


P  = 


E1 


(3.2) 


RT  cavities  typically  have  a  Q  ~  2xl04  and  an  Ea  ~  2  MV/m  compared  to  values  for  SC 
cavities  of  3xl09  and  10  MV/m  respectively.  Some  actual  values  are  shown  in  Table  3.  The 
dissipated  power  is  supplied  to  the  accelerator  by  the  rf  power  source. 


cavity 

/ 
(MHz) 

(MV/m) 

Q 

Z 
(Q/m) 

TESLA  (SC)  [19] 

1300 

20 

>5xl09 

1100 

CEBAF  (SC)  [2] 

1500 

8 

5xl09 

1020 

CDRL  (SO*  [20] 

500 

5.3 

2xl09 

400 

TEUFEL  (RT)  [21] 

1300 

2 

18330 

1714 

Table  3  -  Typical  cavity  parameters  (*  -  parameters  from  a  proposal) 

In  order  to  properly  apply  equation  3.2,  it  is  necessary  to  more  fully  characterize  the 
parameters.  Superconductivity  is  characterized  by  the  BCS  (Bardeen,  Cooper,  and  Schrieffer) 
theory.  For  a  SC  cavity,  the  surface  resistance,  Rs  =  RBCS  +  Rres,  is  given  by  the  relation 
[16][17] 

MT) 

(3.3) 


Rs  =  A  lj  e      B      +  RT 


where  A  (Q  K/m/MHz2)  is  the  BCS  material-dependent  constant,  A  (eV)  is  the  energy  gap  of 
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the  SC  material  used  to  form  Cooper  pairs,  T  (K)  is  the  absolute  temperature  of  the  SC 
structure,  kB  (eV/K)  is  the  Boltzmann  constant,  and  Rres  (Q/m )  is  the  residual  losses  of  the  SC 
structure  after  fabrication  (  Rres  =  1-10  nCl/m  ).  The  quality  factor  of  an  accelerator  cavity  is 
the  ratio  of  the  energy  stored  within  the  cavity  to  the  rf  energy  dissipation  within  the  cavity  per 
rf  cycle.  The  Q  value  is  related  to  the  surface  resistance  through  a  geometry  factor,  G ,  by 

Q=  —  ,  -  =  -*-  +  —  (3.4) 

Rs  Q  QBCS  Qres 

where  QBCs  is  the  contribution  to  the  quality  factor  due  to  losses  in  the  superconductor  and 
Qres  is  the  residual  losses  contribution  to  the  quality  factor.  The  geometry  factor  varies  from 
270  to  290  ohms  [17].  An  average  value  of  G  =  280  Q.  is  assumed.  With  proper  fabrication 
and  conditioning  techniques,  a  Qres,  of  better  than  5xl09  can  be  obtained  [16][17]. 

The  shunt  impedance  is  also  a  function  of  frequency.  In  the  absence  of  other  changes, 
the  shunt  impedance  is  proportional  to  frequency  [22].  However,  as  the  frequency  changes,  the 
subsequent  scaling  of  the  cavity  dimensions  allow  us  to  make  more  favorable,  slight 
adjustments  in  cavity  geometry.   A  curve  fit  to  existing  SC  cavities  is  [17] 

Z  =380(/7500)9        (Q/m).  (3.5) 

The  cavity  geometry  changes  do  not  have  a  very  large  effect  on  the  shunt  impedance.  This  is 
due  to  characteristics  of  SC  cavities  like  the  very  large  apertures  which  can  be  seen  in  Figure 
11.  As  a  result,  the  variations  are  incorporated  into  the  frequency  scaling  of  equation  3.5 
instead  of  the  geometry  factor. 

SC  cavities  are  currently  made  exclusively  of  niobium  (Nb).  Although  there  are  many 
SC  materials  available,  it  is  desirable  to  use  a  material  with  the  highest  superconducting  critical 
temperature,  (Tc ),  and  that  can  support  the  highest  rf  field  strength,  (Hc),  without  forcing  it  out 
of  its  superconducting  state.  It  must  also  have  a  widespread  availability,  relative  ease  of 
fabrication,  high  thermal  conductivity  to  dissipate  heat,  significant  material  phase  stability,  and 
have  stable  superconducting  properties.  Only  lead  (Pb)  and  Nb  currently  meet  these 
requirements  with  Nb  having  the  larger  values  of  Tc  and  Hc.  High  temperature 
superconductors  are  currently  unsuitable  for  these  applications  because  of  large  surface 
resistances  (Rs  £  10~5Q)  and  low  Hc  (the  highest  values  are  1/16  that  of  Nb)  [16]. 

The  BCS  material  constant  for  Nb  is  A  =  O00028/(500)2.  The  static  heat  load  of  a  SC 
cavity  due  to  power  couplers,  etc.  is  typically  about  8  W/m  at  500  MHz  [17].  Therefore,  the 
total  power  dissipated  in  the  cavity  is  PL  =  Pslatlc  +  PBcs  +  pres  >  mat  ls 
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S  E}  Ea2(fl500)1  '  (.000028)?  _1767/r 

PL  =    ,    g      + — +^^ (W/m)  .       (3.6) 

L       V/7500     380(//900)9(2re,  380(280)7 

RT  cavities,  normally  fabricated  from  copper,  have  a  surface  resistance  which  is 
proportional  to  the  square  root  of  the  frequency  and  the  cavity  Q  is  inversely  proportional  to 
the  square  root  of  the  frequency  for  a  fixed  geometry  [22] [23].   Therefore, 

Ea  W 

Pl  =  tV  a7) 

and 

o  -  if  ■  a8) 

where  fc]  and  k2  are  geometry  constants  that  are  inversely  proportional  to  the  square  root  of 
both  the  electrical  resistivity  and  magnetic  permeability  of  the  cavity  material.  The  shunt 
impedance,  then,  is  proportional  to  frequency  and  independent  of  the  material  from  which  the 
cavity  is  fabricated  or  its  surface  finish  [22].  The  constants,  however,  are  heavily  dependent  on 
cavity  geometry.  Unlike  with  SC  structures,  the  geometry  constants  vary  widely  with  the  rf 
frequency  as  adjustments  to  cavity  geometry  are  made  with  the  scaling  of  the  cavity 
dimensions.  This  is  apparent  from  Figure  11  with  the  RT  cavity  cell's  narrow  apertures.  The 
values  of  fcj  for  existing  designs  vary  over  the  range  0.9  <  it,  <  4  (Q/m/MHz1/2) 
[12]  [20]  [21]  [23].   Lower  frequency  cavities  will  have  larger  values  of  k ,. 

The  electron  beam  is  a  series  of  micropulses  that  are  separated  in  space  by  some  integer 
number  of  rf  wavelengths,  \rj .  These  micropulses  can  either  be  injected  into  the  accelerator  in 
macropulses  or  in  CW  mode.  For  a  high-power  FEL,  CW  mode  is  preferred  [13]  because: 

1.  The  peak  power  of  a  micropulse  is  reduced.  In  reducing  the  total  charge  contained 
within  a  micropulse,  it  is  easier  to  overcome  the  disruptive  space  charge  effects  that 
degrade  the  electron  beam  quality.  Furthermore,  the  problem  of  electron  beam  "halo"  is 
reduced.  A  real  electron  beam  has  a  transverse  spatial  distribution.  Although  the  electron 
density  away  from  the  beam  transport  axis  is  a  small  fraction  of  that  on  axis,  the  number 
of  these  electrons  in  a  high  current  pulse  is  significant.  These  electrons  can  strike  the 
magnetic  beam  transport  optics  and  deflect  so  that  they  excite  rf  modes  and  cause 
transport  instabilities  for  the  entire  beam.  Similarly,  the  activation  of  the  optics  would 
create  a  significant  radiation  hazard  over  a  very  short  time  [12]  [14]. 


33 


2.  The  peak  power  of  the  optical  pulses  developed  in  the  optical  resonator  is  reduced 
extending  the  lifetime  of  the  mirrors  and  various  optical  components.  Furthermore,  the 
more  constant  loading  of  the  optics  reduces  the  need  for  stabilization  measures  due  to 
transient  deformations. 

3.  The  overall  power  efficiency  of  the  accelerator  is  increased.  When  electron  pulses  are 
accelerated  in  macropulses,  the  rf  power  is  applied  to  the  accelerator  in  pulses  (pulsed 
power  operation).  A  finite  "filling  time"  is  required  to  build  up  the  energy  storage  within 
the  accelerator  cavity  prior  to  the  injection  of  the  next  macropulse.  There  is  significant 
energy  lost  in  the  accelerator  during  both  this  filling  time  and  the  decay  time  at  the  end 
of  the  macropulse  where  stored  energy  within  the  cavity  is  dissipated. 

4.  Optimization  of  an  energy  recovery  system  is  possible.  Pulsed  power  systems  have 
significant  time  variations  of  the  electron  beam  parameters  over  the  duration  of  the 
macropulse.  This  affects  the  optical  evolution  in  the  resonator.  In  CW  mode,  these 
variations  occur  over  only  a  small  fraction  of  the  electron  beam's  duration  and  are  not 
repetitive. 

5.  Finally,  the  control  systems  are  easier  to  realize.  They  can  be  simpler,  slower,  and  much 
less  expensive  than  for  an  equivalent  pulsed  power  machine.  This  is  an  extremely 
important  factor  because  the  feedback  and  control  requirements  for  a  high-average  power 
energy  recovery  system  are  formidable.  Studies  into  the  control  system  requirements  for 
these  systems  are  just  beginning  [10]  [12]. 

The  final  consideration  in  determining  the  power  requirements  is  in  the  choice  of  rf 
frequency.  In  general,  it  is  preferred  to  have  as  low  a  frequency  as  possible  without  the  size 
and  capital  cost  of  the  accelerator  structure  becoming  excessive.  Lower  frequency  cavities  are 
preferred  because  [13] [22]  [24]: 

1 .  The  highest  average  power  rf  sources  are  at  500  MHz  and  below. 

2.  For  CW  applications,  which  would  become  necessary  for  high-average  power  electron 
beams,  cooling  considerations  for  RT  cavities  favor  lower  frequencies.  Although  the  heat 
losses  are  larger  at  lower  frequencies,  the  larger  cell  surface  area  provides  for  more 
effective  heat  transfer. 

3.  The  larger  surface  area  per  cell  allows  for  more  effective  and  simpler  cooling  of  the  cell 
walls. 

4.  Lower  frequency  cavities  have  the  potential  to  produce  higher  quality  electron  beams  for 
three  reasons.    First,  for  a  given  spatial  length  of  each  electron  micropulse,  it  will  span  a 
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smaller  A(J>  phase  of  the  rf  accelerating  fields  which  leads  to  more  uniform  acceleration 
properties.  Second,  the  transverse  dimensions  (apertures)  of  the  cavity  are  larger  so  the 
"transverse  wake-field  effects,"  or  beam  emittance  induced  in  traveling  through  an 
aperture  due  to  the  image  charges  formed  on  the  surface,  is  reduced.  Third,  the  cells  are 
physically  longer  since  the  length  of  a  cell  is  approximately  Xrf  12.  Therefore,  for  a  given 
acceleration  gradient,  the  beam  will  pass  through  fewer  apertures  to  reach  a  given  energy. 

5.  The  volume  of  the  cavity  is  proportional  to  X^  so  a  reduced  frequency  leads  to  an 
increased  cavity  volume  and  higher  stored  energies.  The  electron  beam  traveling  through 
the  cavity  creates  higher  order  modes  (HOM)  by  virtue  of  its  charge.  In  higher  energy 
density  cavities,  these  HOM  will  disturb  the  cavity  fields  less  and  reduce  the  threshold  of 
beam-breakup  instabilities  that  can  result. 

6.  Individual  cavities  can  be  more  easily  electromagnetically  uncoupled  which  allows  for 
better  stabilization  of  the  fields  within. 

The  disadvantages  of  lower  frequency  operation  are  the  increased  size  of  the  structure, 
and  the  fact  that  acceleration  gradients  must  be  lower.  They  are  also  more  difficult  to 
manufacture  because  of  the  increased  size. 

Figure  12  is  a  plot  of  equation  3.6  given  values  of  Ea  =  10  MV/m,  W  =  E  x  La  =96 
MeV,  Qres  =  8xl09,  and  T  =  4.2  K.  The  minimum  occurs  at  a  frequency  of  300  MHz.  If  Qres 
or  T  decreases  slightly,  then  this  minimum  will  shift  to  higher  frequencies.  Therefore,  a 
frequency  of  400  MHz  is  a  suitable  design  point  for  the  SC  accelerator.  The  wall  losses  for  the 
SC  structures  will  be  approximately  PL  =  1.2  kW.  This  value  is  insignificant  to  the  overall 
power  efficiency  of  the  system.  However,  this  value  will  determine  the  requirements  for  the 
liquid  helium  refrigeration  support  system. 

The  RT  accelerator  operating  frequency  is  not  as  easy  to  determine.  There  are  three 
criteria  which  determine  the  possible  operating  frequency  range.  First,  the  maximum  field 
strength,  £max,  allowed  within  a  cavity  before  sparking  occurs  is  £max  =  \.5EKPIDM  [23] 
where  DM  is  an  appropriate  design  margin  and  E^  (MV/m)  is  the  Kilpatrick  sparking  criteria 
field  strength  determined  from  the  trancendental  relation 

fmlM  Ej>  €-***».  (3.9) 

A  design  margin  of  four  is  generally  used. 

Second,  a  minimum  field  strength  within  the  cavity  is  necessary  to  ensure  adequate 
energy  storage  so  that  the  energy  removed  in  accelerating  an  electron  pulse  is  a  small  fraction 
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SC  -  Wall  Losses 
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Figure  12  -  SC  accelerator  wall  losses 

of  the  total  cavity  stored  energy.   This  minimum  field  strength,  £min  is  given  by  [23] 

_  k]  2nf  M  lavg 
min      k2      (AUIU) 


(3.10) 


where  AU IU  is  the  fraction  of  energy  removed  from  the  cavity  by  an  accelerating  electron 
pulse  and  M  represents  the  pulse  repetition  rate  in  integer  numbers  of  rf  cycles.  A  value  of 
M  =  1 2  means  that  an  electron  pulse  is  accelerated  every  1 2th  rf  cycle. 

The  final  criteria  needed  to  determine  the  operating  frequency  for  RT  structures  is  the 
maximum  allowed  wall  losses,  Pmax.  Equation  3.7  establishes  the  relationship  between  Ea  and 
/  for  a  fixed  P L . 

The  relationship  between  these  three  criteria  is  shown  in  Figure  13(a)  with  values 
A,  =  3.5  MWmlMHz"2,  k2  =  8xl05  MHzm,  M  =  12,  Iavf,  =  500  mA,  AUIU  =  0.02,  and 
Pmax  =  4  MW.  The  Ea  .    curve  can  be  lowered  only  by  changing  M    given  the  design 

constraints  of  our  problem.  Current  accelerators  operate  with  values  of  12  <  M  <  20.  For  a 
given  value  of  M  and  AUIU ,  the  wall  losses  are  less  at  lower  frequencies  because  of  the 
reduction  in  field  strength.  A  value  of  M   =  12  was  chosen  to  be  consistent  with  current 
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Acceleration  Gradient  Limitations 
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Figure  13  -  RT  accelerator  (a)  field  strength  requirements,  and  (b)  wall  losses 

technologies.    As  M   is  reduced,  the  peak  currents  are  reduced  as  is  the  excitation  of  HOM 
within  the  accelerator  cavities.   However,  less  time  is  available  to  damp  these  modes  out. 

Within  the  design  constraints  of  the  FEL,  a  maximum  rf  frequency  is  desired  to  limit  the 
wall  losses  in  the  RT  structure.  A  maximum  acceleration  gradient  is  also  desired  in  order  to 
reduce  the  active  length  of  the  accelerator  structure.  However,  this  is  limited  by  the  resultant 
increased  power  losses.  Therefore,  an  operating  frequency  of  200  MHz  is  chosen  for  the  RT 
accelerator.    From  Figure  13(b),  the  wall  losses  are  PL  =  2.9  MW  with  an  Ea  =  1.5  MV/m. 

This  analysis  does  not  address  the  suitability  of  a  structure  with  these  parameters  for  use 
in  an  FEL.  Extensive  modeling  and  numerical  simulations  would  be  required  to  determine  the 
beam  quality  and  beam  transport  characteristics  that  could  be  achieved  with  such  a  structure. 
However,  these  parameters  have  been  validated  by  a  comparison  to  an  actual  180  MHz  cavity 
design  [12]. 

It  is  prudent  to  note  that  the  trend  for  future  compact  FEL  accelerators  is  to  include 
recirculation  with  energy  recovery  where  a  beam  is  bent  around  to  pass  through  the  accelerator 
two  or  more  times  to  reach  its  final  energy  [2][  13][25).  This  allows  accelerators  to  be  more 
compact  in  that  the  active  length  of  the  accelerator  is  smaller  by  a  factor  approximately  equal 
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to  the  number  of  passes.    The  power  estimates  for  these  specific  architectures  will  be  nearly 
identical  to  the  analysis  presented  here.   However,  the  control  system  issues  are  more  complex. 

C.  ELECTRON  INJECTOR 

The  purpose  of  the  electron  injector  is  to  produce  high  density  pulses  of  electron  current 
and  inject  them  into  the  accelerator  at  the  proper  phase  relation  to  the  accelerating  rf  fields. 
The  injector  consists  of  two  major  components:  the  electron  gun  and  the  buncher.  The  buncher 
systems  contain  rf  cavities  similar  to  accelerator  cavities.  These  cavities  are  designed  to 
produce  bunched  electron  pulses  of  low  emittance.  One  of  the  most  promising  electron  gun 
technologies  for  the  production  of  high  current,  low  emittance  electron  beams  center  on  the  dc 
photocathode  [12]  [14]. 

Figure  14  shows  the  functional  component  interconnections  of  an  injector  driven  by  a  dc 
photocathode  electron  gun.  When  illuminated  by  laser  light,  electrons  are  liberated  at  the 
photocathode  through  the  photoelectric  effect.  A  dc  high-voltage  rapidly  accelerates  the 
electrons  away  from  the  photocathode  and  directs  them  into  the  buncher.  The  dc  voltage 
should  be  large  to  prevent  the  Coulombic  space  charge  effects  from  inducing  excessive  initial 
emittance.  However,  since  the  electron  velocity  cannot  exceed  the  speed  of  light,  the  voltage 
must  be  low  enough  to  limit  the  electron  velocity  sufficiently  below  the  speed  of  light  when 
the  electrons  enter  the  buncher  so  it  can  effectively  bunch  the  micropulses  before  they  enter  the 
accelerator  [22].  A  suitable  value  for  the  dc  voltage  is  500  kV.  The  electrons  have  been 
accelerated  to  a  kinetic  energy  of  about  500  keV  (  velocity,  v  «  0.85c)  when  they  enter  the 
buncher. 

In  describing  the  electron  gun  power  requirements,  it  is  more  convenient  to  start  at  the 
high  voltage  dc  electrodes  and  work  back  towards  the  drive  laser.  The  high  voltage  has  a 
practical  limit  due  to  dielectric  breakdown  effects  and  the  emittance  limiting  effects  previously 
mentioned.  As  buncher  and  high  voltage  power  supply  technology  has  improved,  however, 
this  voltage  limit  has  also  risen.  Currently,  and  for  the  near  future,  these  devices  can  be 
expected  to  operate  at  about  500  kV  [2] [14].  Therefore,  for  a  500  mA  average  electron 
current,  250  kW  of  high  voltage  dc  electrical  power  must  be  supplied  by  the  dc  electric  field  to 
the  electron  current. 

To  generate  the  500mA  of  average  current,  we  use  a  photocathode  which  emits  electrons 
via  the  photoelectric  effect.  It  is  assumed  that  our  accelerator  will  be  operating  in  the  CW 
mode  so  that  the  peak  currents  in  each  micropulse  is  minimized.  However,  to  minimize  the 
beam  emittance,  the  electron  micropulses  should  have  a  spatial  length  corresponding  to  only  a 

38 


+ 


Nd  Laser 


Pulse 

Compression 


Optical 
Amplifier 


Frequency 
Doubling 


Photo- 
Cathode 


(a)  typical  dc  photocalhode  electron  gun 


Electron 
Gun 


Buncher 

C)C) 


Accelerator 


(b)  electron  injector 


Figure  14  -  typical  electron  injection  system 


few  degrees  of  rf  phase.  Reference  [22]  defines  a  good  injector  as  supplying  a  phase  spread  of 
A<J>  £  10°  and  a  very  good  injector  as  having  a  A<J>  <  5°.    Low  emittance  electron  beams  are 

being  produced  with  micropulse  phase  spreads  of  approximately  A<}>  =  5-6°  [2][14][26].  For  a 
CW  400  MHz  system,  a  A<|>  =  6°  corresponds  to  a  pulse  length  of  42  ps  with  an  average 
micropulse  current  of  30  A  per  rf  cycle.  Therefore,  we  require  an  average  electron  charge  of 
1.25  nC  per  rf  cycle  or  1.25xAf  nC  per  micropulse.  At  a  rf  frequency  of  200  MHz,  we 
require  an  average  electron  charge  of  2.5xA/  nC  per  micropulse.  Electron  pulses  of  3-4  nC, 
peak  currents  of  several  hundred  amps,  and  average  currents  of  tens  of  mA  have  been 
effectively  produced  [26][2][14]  with  rapid  advancements  being  made.  The  expectation  of  a 
500  mA  average  current  in  the  future  is  reasonable. 

The  photoelectric  effect  is  described  by 


El„  =  h  v  +  W 


-i, 


r, 


(3.11) 
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where  Eke  is  the  kinetic  energy  of  the  liberated  electron,  hv  is  the  energy  of  the  absorbed 
photon,  and  W^  is  the  work  function  of  the  photoemissive  material.  The  values  of  the  work 
functions  for  the  different  elements  range  from  3.14  eV  (cesium)  to  about  6  eV  [27].  This 
minimum  photon  energy  required  to  generate  a  free  electron  means  that  photocathodes  are 
insensitive  to  infrared  radiation.  The  amount  of  charge  liberated  by  a  pulse  of  total  energy  Ep 
with  photons  of  energy  h  v  is 


q(nC)=  10 


Ep(\ti)QE{%) 


(3.12) 


/jv(eV) 

where  QE  is  the  effective  quantum  efficiency  which  is  the  ratio  of  the  number  of  liberated 
electrons  to  the  number  of  incident  photons.  Table  4  lists  the  effective  quantum  efficiency  for 
a  number  of  photocathodes  under  development. 


photocathode 

QE 
(%) 

"•incident 

(nm) 

Comments 

CsK  2Sb  [26] 

0.1 

527 

1-3  nC 

Cu  [UCLA] 

.01 

266 

>  1  nC 

Ag  [29] 

.03 

248 

Cu+10%BaO  [Aganofov] 

.17 

266 

LaB6[31][32] 

.06 

349 

T=  1500  K 

GaAs  [14]  [2] 

18 

527 

CsTe  [33]  [34]  [35] 

13 

263 

long  lifetime 

CsTe  [34]  [35] 

19 

251 

Table  4  -  Photocathode  quantum  efficiencies 

It  is  reasonable  to  expect  that  a  photocathode  will  be  developed  suitable  for  a  shipboard  FEL 
with  a  quantum  efficiency  of  QE  =  20%.  However,  a  QE  =  5%  would  be  a  more 
conservative  expectation,  and  would  account  for  QE  variations  over  the  lifetime  of  the 
photocathode.  From  3.12  we  can  show  that  the  required  laser  input  power  to  the  photocathode 
to  generate  a  required  current  is 
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P  =     /?v(eV)     /(mA)  (W).  (3.13) 

10Q£(%) 

For  a  frequency  doubled  Nd:YAG  laser  (A  =  532  nm)  and  a  (2£  ~  5%,  an  average  current  of 
500  mA  would  require  an  input  power  of  1 1 .7  W  incident  on  the  photocathode. 

The  frequency  doubler  is  a  non-linear  electrooptic  crystal  in  which  second  harmonic 
generation  occurs.  If  we  align  the  crystal's  crystallographic  axes  such  that  the  index  of 
refraction  for  both  the  fundamental  and  second  harmonic  frequencies  are  equal,  the  crystal  is 
"phase-matched".  Under  the  phase-matched  condition,  two  photons  of  the  fundamental 
essentially  combine  to  form  a  new  photon  of  twice  the  frequency,  that  is 

Tito  +  Tito  =  7i(2(d)  (3.14) 

The  efficiency  of  frequency  doubling,  ranges  from  20-65%  depending  on  the  input  power  level. 
Present  photocathode  electron  guns  operate  with  a  frequency  doubling  efficiency,  x\fd ,  of  about 
50%  [26]  [14],  but  that  can  vary  depending  on  the  actual  incident  power  intensity.  For  linac 
applications,  the  length  of  the  electro-optic  crystals  are  limited  by  group  velocity  mismatch 
which  tends  to  increase  the  pulse  length  and  limit  the  frequency  doubling  efficiency.  Using 
type  1  phase-matching,  however,  efficiencies  of  more  than  60%  can  be  obtained  [36].  Due  to 
our  high-power  requirements,  a  frequency  conversion  efficiency  of  J\fd  =  55%  is  assumed. 

An  optical  amplifier  is  a  pumped  laser  medium  which  amplifies  the  input  signal  through 
stimulated  emission.    An  optical  amplifier  is  described  by  the  expression  [37] 


*(0) 


ln(Gsp)+^-[Gsp-l]  =  V0d  (3.15) 


where  (30  is  the  unsaturated  gain  coefficient,  d  is  the  single  pass  length  of  the  amplifier 
medium,  Gsp  is  the  overall  single  pass  gain,  <t>s  is  the  saturation  photon  flux  density,  and  4>(0) 
is  the  input  photon  flux  density.  <t>s  is  a  function  of  only  the  amplifier  medium  while  f30  is  a 
function  of  the  amplifier  medium  and  is  directly  proportional  to  the  population  inversion  per 
unit  volume.   Therefore,  f30  will  determine  the  power  requirements  for  the  amplifier. 

To  maximize  the  power  efficiency  of  the  amplifier,  the  amplifying  medium  must  be 
efficiently  pumped  and  multiple  passes  must  be  made  by  the  light  wave  through  the  amplifying 
medium.  To  efficiently  pump  the  medium,  laser  diode  pumping  provides  the  most  efficient 
method  because  the  pump  energy  is  concentrated  at  the  pump  bands  of  the  Nd  lasing  medium. 
The  width  of  the  diode  pump  spectrums  are  on  the  order  of  5  nm  FWHM  [38].    Furthermore, 
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diode  lasers,  driven  by  a  low  dc  voltage,  can  achieve  a  quantum  efficiency  of  near  unity 
[39]  [40].  One  particular  diode-pumped  optical  amplifier  design  utilizing  a  confocal,  parabolic 
mirror  cavity  with  an  internal  Brewster  plate  to  increase  the  path  length  of  the  beam  within  the 
crystal,  has  achieved  small  signal  gains  of  38  dB  with  a  pump  to  optical  power  efficiency  of 
40%  for  the  Nd:YLF  medium  and  25%  in  the  Nd:YAG  medium  operating  in  the  CW  mode 
[41].  Amplified  spontaneous  emissions  are  only  about  1%  of  the  total  output  power.  The 
Nd.YAG  amplifier,  however,  was  able  to  support  much  more  pump  power  than  the  YLF 
amplifier  which  fractured  at  lower  pump  powers.  For  the  estimation  of  amplifier  power 
requirements,  it  is  assumed  that  an  optical  amplifier  for  my  FEL  photoinjector  has  a  pump 
power  amplification  efficiency,  i\amp ,  of  25%. 

Nd:YLF  and  Nd:YAG  mode-locked  CW  lasers  are  commonly  used  as  photoinjector  drive 
lasers  because  they  are  relatively  rugged,  produce  short  pulses  (about  75  ps  for  YAG  [39]  and 
50  ps  for  YLF  [42]),  have  large  optical  gain  per  unit  distance,  and  can  be  easily  amplified. 
The  large  optical  gain  means  they  have  a  low  threshold  power  required  for  CW  operation. 
With  laser  diode  pumping,  we  can  achieve  much  higher  pumping  efficiencies.  The  resulting 
inefficiencies  are  predominantly  due  to  incoupling  losses,  the  fraction  of  energy  absorbed  by 
the  lasing  medium,  and  the  spontaneous  emission  losses.  Optical-to-optical  efficiencies  of  the 
diode-pumped  laser  have  been  demonstrated  at  25%  for  Nd  lasers  in  CW  mode  [43].  Diode 
lasers  delivering  60  W  for  a  200  usee  duration  are  commercially  available.  Advances  in  the 
development  of  higher-powered  diode  lasers  and  improvements  in  the  coupling  optics  are 
expected  to  lead  to  pulse  durations  below  10  nsec  in  the  near  future  [43].  This  is  sufficient  to 
support  the  42  ps  pulse  lengths  needed  to  support  this  FEL. 

Current  Nd  drive  lasers  used  in  photocathode  electron  guns  have  self  phase  modulation 
(SPM)  pulse  compression  systems  with  an  incoupling  loss  of  25%  [42].  These  guns  drive 
accelerators  operating  at  rf  frequencies  of  several  GHz  where  narrow  electron  pulses  are 
required.  The  pulse  compression  systems  do  not  require  input  biasing  or  pumping  power  but 
do  limit  the  peak  output  power  from  the  laser  to  about  2  W.  Pulse  compression  systems, 
however,  should  not  be  required  for  the  low  accelerator  frequencies  under  consideration  for  a 
high-power  FEL.  If  we  assume  a  laser  output  power  of  2  W  at  a  wall  plug  efficiency  of  25%, 
only  a  few  watts  of  power  are  needed.  Therefore,  the  drive  laser  power  and  pulse  compression 
losses  can  be  neglected  from  the  prime  power  analysis. 

The  power  requirements  of  the  electron  gun  now  can  be  summarized  in  terms  of  Figure 
15.  The  total  required  low  voltage  dc  power,  FWf,  is  given  by 
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Figure  15  -  Power  flow  in  the  dc  photocathode 
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and  the  total  high  voltage  dc  required,  P^dc,  is  the  average  current  times  the  voltage  applied. 
For  a  500  mA  average  current  and  500  kV  acceleration,  we  will  require  PWc  =  170  W  and 
PHdc  =  250  kW. 

The  buncher  cavities  are  RT  cavities  whose  geometry  and  fields  are  optimized  to  provide 
bunching  of  the  electrons  into  micropulses  with  low  emittance  and  energy  spread.  SC  cavities 
and  their  associated  geometries  are  not  efficient  for  tins  purpose  so  RT  cavities  are  used 
exclusively.  The  electrons  liberated  from  the  photocathode  are  bunched  longitudinally  and 
focused  to  a  waist  as  they  are  injected  into  the  accelerator.  Solenoidal  magnetic  fields  are  also 
used  to  assist  the  bunching  process. 

For  the  defined  architecture  and  parameters,  the  buncher  cavities  must  provide  3  MeV 
acceleration  to  the  electrons  during  this  bunching  process.  Since  the  electrons  are  entering  the 
buncher  at  a  velocity  of  v  =  0.85c,  the  fields  must  be  rather  strong  to  bunch  the  electrons.  To 
determine  the  wall  losses  in  the  buncher,  PLi„r  the  same  model  used  for  estimating  the 
accelerator  wall  losses  is  employed.  To  achieve  3  MeV  acceleration,  three  cavity  cells  at  a 
frequency  of  200  MHz  and  Ea  =  1.33  MV/m  is  assumed.  To  account  for  added  losses  due  to 
the  required  geometry  changes,  the  shunt  resistance  constant  is  assumed  to  be  A,  =  1.9 
(Mft/m/MHz,/2).  Therefore,  we  can  estimate  PLinj  a  150  kW.  At  a  frequency  of  400  MHz, 
only  100  kW  of  power  would  be  dissipated. 

Although  electromagnets  are  often  used  to  produce  the  solenoidal  fields,  large  permanent 
magnets  with  small  tuning  electromagnets  would  be  more  appropriate  for  high-current 
applications.   Therefore,  they  will  not  add  significantly  to  the  power  requirements. 
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D.  RF  SYSTEMS 

Electron  accelerators  have  historically  operated  at  rf  frequencies  of  more  than  1  GHz,  low 
average  power,  and  high  peak  powers  in  a  pulse  power  electron  beam  format.  Klystrons  have 
been  the  only  rf  sources  suitable  for  these  operating  conditions.  For  rf  frequencies  of  200-500 
MHz,  however,  additional  rf  sources  may  be  practical. 

Table  5  is  a  listing  of  several  of  the  most  powerful  CW  klystrons  currently  on  the 
market. 


TH2089 

YK1350 

YK1303 

K3513 

Units 

Reference 

[44] 

[44] 

[44] 

[45] 

frequency 

352 

352 

508 

352 

MHz 

input 
voltage 

87 

90 

90 

100 

kV 

input 
current 

17.1 

16.8 

18.2 

20 

A 

Power 

1 

1 

1 

1.3 

MW 

efficiency 

68 

66 

61 

65 

% 

length 

4.75 

3.75 

3.75 

4.2 

m 

Table  5  -  Commercially  available  super-power  CW  klystrons 

Klystrons  have  evolved  over  several  decades  so  revolutionary  advances  in  their  power 
capabilities  within  the  next  20  years  are  not  anticipated.  Some  of  these  tubes  are  designed  to 
operate  in  a  vertical  orientation  while  others  are  designed  to  operate  in  a  horizontal  orientation. 
High-power  klystrons  are  available  in  either  configuration.  The  output  power  and  power 
conversion  efficiency  of  a  klystron  for  a  shipboard  FEL  can  be  expected  to  have  an  output 
power  of  about  1.5  MW  at  an  efficiency  of  65%. 

Some  alternative  rf  sources  for  accelerator  applications  are  listed  in  Table  6.  The  high- 
power  entries  in  this  table  are  projections  and  have  yet  to  be  developed.  However,  these 
alternate  technologies  promise  rf  sources  that  operate  at  lower  voltages,  higher  efficiency,  and 
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MBK 

MBK* 

klystrode 

klystrode* 

IOT 

Units 

Reference 

[46] 

[46] 

[47] 

[48] 

[49] 

frequency 

425 

500 

267 

476 

476 

MHz 

input 
voltage 

17 

39 

65 

25 

66 

kV 

input 
current 

11 

41 

5.5 

70 

9.3 

A 

Power 

0.1 

1 

0.25 

1.3 

0.4 

MW 

efficiency 

55 

64 

70 

75 

>65 

% 

length 

1.35 

2.5 

NA 

NA 

2.1 

m 

400 

800 

NA 

NA 

500 

kg 

Table  6  -  Alternative  CW  rf  power  sources  (*  -  denotes  projection  and  is  not 
commercially  available) 


have  a  reduced  component  weight  and  volume.  It  is  expected  that  the  multiple  beam  klystron 
(MBK)  efficiency  can  be  increased  to  >80%  with  further  development  [24].  The  klystrode  and 
inductive  output  tube  (IOT)  are  essentially  the  same  device  with  klystrode  being  Varian 
Associates  trademark  for  this  device.  It  has  the  disadvantage  of  a  smaller  gain.  Therefore,  a 
larger  rf  drive  is  required  for  the  device  (several  kW  for  a  megawatt  device).  This  will  add  to 
the  complexity  of  the  overall  rf  system. 

E.  SUPPORT  SYSTEMS 

A  wide  variety  of  support  systems  will  be  required  for  a  shipboard  FEL.  Feedback 
control  systems  will  be  required  to  maintain  rf  power  source  phase  and  amplitude  stability.  The 
2.5  MW  of  heat  energy  dissipated  in  the  beam  dump  and  the  3  MW  of  heat  energy  dissipated 
in  the  RT  structure  as  wall  losses  will  need  an  appropriate  cooling  system.  System  monitoring 
and  fire-control  interconnections  will  also  place  a  load  on  the  ship's  electrical  system.  These 
systems,  however,  can  easily  be  adapted  to  existing  shipboard  systems.  Current  steam 
propulsion  plants  dump  heat  energy  well  in  excess  of  50  MW  at  a  cost  of  tens  of  kW  of 
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seawater  pump  power.  Therefore,  a  reasonable  estimate  is  that  these  systems  will  add  an 
additional  100  kW  to  the  ship's  electrical  load. 

There  are  three  support  systems,  however,  that  need  to  be  addressed  separately.  First,  a 
SC  structure  will  require  a  liquid  helium  refrigeration  system.  The  power  requirements  for  this 
system  can  be  estimated  by  [17] 

T„  -THe      i 
Pr    =  Pl  ~V^  ~r  (3.17) 

THe  T\R 

where  PR  is  ship's  service  AC  power  consumed  by  the  refrigeration  system,  Ta  is  the  ambient 
(room)  temperature  (Ta  ~  300/O,  THe  is  the  liquid  helium  boiling  point  temperature 
(THe  ~4.2K),  and  T]R  is  the  refrigeration  efficiency.  The  second  term  in  3.17  represents  the 
Carnot  thermodynamic  efficiency  of  the  system  and  the  third  term  is  the  refrigeration  system 
(compressor,  etc.)  efficiency.  Typical  values  of  r\R  are  approximately  25-35  percent  [14][16]. 
For  a  30%  efficient  refrigeration  system  removing  a  PL  =  1.2  kW,  our  shipboard  FEL  will 
require  PR  ~  (235xPL )  =  250  kW  of  refrigeration  power  if  the  laser  were  to  operate 
continuously.  This  is  not  only  a  significant  amount  of  power  but  the  size  of  the  refrigeration 
plant  would  be  large  [2]  [14]. 

The  refrigeration  plant  for  this  shipboard  FEL,  however,  can  be  considerably  smaller 
since  it  will  not  operate  continuously.  The  rate  of  liquid  helium  boiled  off  during  FEL 
operation  is 

V  =  -^—  PL  (3.18) 

P  <7v 

where  V  is  the  rate  of  helium  boiled  of  during  laser  operation  (gal/sec),  A  is  the  atomic  weight 
of  helium  (4  g/mole),  p  is  the  density  of  helium  at  THe  (7.62  lb/ft3),  and  qv  is  the  latent  heat 
of  vaporization  of  helium  (84  J/mole).  The  operational  requirements  specify  that  the  FEL  must 
operate  for  150  seconds  over  a  20  minute  interval.  Conservatively  assuming  that  the  laser 
operates  continuously  for  150  seconds  with  a  wall  loss  of  PL  =  1.2  kW,  18.5  gal  of  liquid 
helium  is  boiled  off.  This  is  a  relatively  small  amount  of  required  stored  liquid  helium. 
Therefore,  the  refrigeration  system  can  be  designed  to  recompress  and  refrigerate  this  boil-off 
over  a  longer  period. 

The  refrigeration  system  sizing  criteria,  therefore,  is  based  on  two  factors.  First,  the 
refrigeration  system  must  be  able  to  remove  the  latent  heat  of  vaporization  from  the  helium 
boiled  off  during  laser  operation  within  a  specified  time  period.  Arbitrarily  choosing  this  time 
period        as        2        hours,        then        the        refrigeration        power        required        is 
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PR  =  250  kW  (100  sec/2  hrs)  =  3.5  kW.  Most  importantly,  the  refrigeration  system  size  and 
weight  has  been  reduced  significantly.  A  reserve  storage  capacity  of  liquid  helium  can  reduce 
this  requirement  even  further.  The  second  factor  is  the  time  required  to  reduce  the  temperature 
of  the  accelerator  from  room  temperature  to  operating  temperature  following  an  import 
maintenance  period.  To  perform  this  cooldown  within  48  hours  would  require  a  refrigeration 
system  that  would  be  excessively  large  for  the  shipboard  environment  [14].  A  shore  facility 
could  perform  this  cooldown  prior  to  a  warship  getting  underway  and  eliminate  the  need  for 
this  sizing  criteria. 

The  second  support  system  that  needs  to  be  addressed  is  the  network  of  electron  beam 
bending  magnets  which  control  the  beam  trajectory.  Large  permanent  magnets  would  also  be 
appropriate  for  this  purpose.  Although  they  are  much  more  expensive,  they  require  no  external 
electrical  power  source.  If  electromagnets  were  used  and  a  magnet  failed,  the  beam  would 
impact  the  wall  of  the  beam  tube  creating  a  large  radiation  hazard  for  adjacent  personnel  and 
causing  a  shutdown  of  the  FEL. 

The  final  support  system  deals  with  the  optical  resonator.  The  resonator  will  have 
megawatts  of  power  circulating  between  the  mirrors.  Special  mirror  configurations  have  been 
designed  for  high-power  applications  [2]  to  reduce  the  power  intensity  on  the  mirrors  and  the 
onset  of  damage.  The  optics  will  need,  however,  an  active  stabilization  system  to  maintain 
alignment  during  ship  vibrations,  flexing,  and  torsional  movements.  Furthermore,  neglecting 
outcoupling  of  the  laser  light,  the  mirrors  will  also  absorb  some  energy  and  need  a  cooling 
system.  Mirrors  have  been  developed  with  a  reflectance  of  99.975%  at  X  =  .87  p.m  [50].  In 
general,  mirror  reflectivity  improves  as  the  optical  wavelength  moves  into  the  infrared. 
Assuming  a  reflectance  of  99.9%  with  a  resonator  stored  power  of  5  MW,  5  kW  of  energy  is 
absorbed  in  the  mirrors.  Suitable  stabilization  and  cooling  systems  have  already  been  developed 
[24][51].  The  power  requirements  can  be  included  with  the  overall  support  system  power 
allocation,  Pss,  of  100  kW. 

F.  SUMMARY  OF  POWER  REQUIREMENTS 

The  power  requirements  for  the  FEL  architecture  of  Figure  10  under  the  design 
constraints  of  a  naval  application  are  summarized  in  Table  7.  Since  about  95%  of  the  electron 
beam  power  is  recovered  in  this  architecture,  the  listed  power  requirements  are  nearly 
independent  of  the  output  power  level.  For  example,  the  accelerator  wall  losses  are  solely  a 
function  of  the  stored  cavity  field  strength.  Therefore,  the  wall-plug  efficiency  of  the  FEL 
scales  as  the  output  power.    The  rf  source  efficiency  is  based  on  commercially  available 
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parameter 

SC  value 

RT  value 

units 

/ 

400 

200 

MHz 

Wb 

100 

100 

MeV 

'avg 

500 

500 

mA 

Ea 

10 

1.5 

MV/m 

Z 

311 

990 

Q/m 

vJres 

5xl09 

Qo 

3xl09 

5xl04 

Pbdc 

350 

350 

W 

PHdc    @  500  kV 

250 

250 

kW 

P 

2 

2 

MW 

Punj 

150 

150 

kW 

p 

1  beam 

50 

50 

MW 

n0 

2 

2 

% 

'oul 

1 

1 

MW 

p 

'  ret 

49 

49 

MW 

'xrec 

95 

95 

% 

p 

'  dump 

2.45 

2.45 

MW 

Prec     @   T\rec    =    95 

46.6 

46.6 

MW 

Pl 

.0012 

2.9 

MW 

pr  @r\R=  -3 

275 

kW 

Prf    @  t\rf  =  -65 

6.6 

9.2 

MW 

p 

100 

100 

kW 

wall-plug  efficiency 

17.4 

10.4 

% 

Table  7  -  Summary  of  the  power  requirements  for  a  shipboard  FEL 

klystrons.  Alternative  rf  sources  promise  an  increased  efficiency  but  require  further 
development  before  they  can  produce  high-average  power.  If  extrapolations  are  made  to  higher 
power  for  these  devices,  then  the  FEL  total  power  requirements  can  be  reduced  by  about  10%. 

For  the  RT  accelerator,  it  was  necessary  to  define  the  pulse  repetition  interval  represented 
by  the  parameter  M  to  determine  the  accelerator  operating  frequency  and  wall  losses.  From 
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equation  3.10  and  Figure  13  it  can  be  seen  that  reducing  the  value  of  M  can  reduce  the  RT 
cavity  wall  losses.  However,  the  optical  resonator  must  have  dimensions  set  so  that  when  an 
electron  pulse  passes  through  the  undulator,  it  overlaps  with  an  optical  pulse  and  can  transfer 
energy  to  the  optical  field.  Therefore,  the  mirrors  must  be  separated  by  a  distance 
L  =  nAfc/2/,  where  n  is  an  integer  and  equals  the  number  of  optical  pulses  stored  within  the 
resonator  cavity.  For  the  values  in  Table  7  and  assuming  M  =  12,  L  =  9n  meters  for  the  RT 
accelerator  parameters.  If  the  resonator  cavity  was  18  meters  long,  then  there  would  be  two 
optical  pulses  circulating  within  the  resonator. 

For  the  SC  accelerator  based  FEL,  M  is  not  an  important  parameter  in  determining  power 
losses  or  the  rf  frequency  because  of  the  high  cavity  Q  and  large  accelerating  field  strengths. 
However,  the  pulse  repetition  rate  will  affect  the  length  of  the  resonator  in  the  same  way  as  for 
the  RT  accelerator  based  FEL.  For  the  SC  parameters  defined  in  Table  7,  and  assuming 
M  =  12,  the  resonator  length  is  L  =  4.5n  meters. 
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IV.  PRIME  POWER  DISTRIBUTION 

In  the  previous  chapter,  the  power  requirements  for  a  megawatt  FEL  laser  weapon  for 
shipboard  applications  were  outlined.  The  two  components  that  burden  the  prime  power 
distribution  system  are  the  dc  photocathode  high  voltage  and  the  rf  sources  for  the  electron  rf 
accelerator.  It  was  determined  that  the  dc  photocathode  would  require  250  kW  at  500  kV  dc 
while  the  rf  sources  would  require  6-10  MW  at  80-100  kV  dc.  This  amount  of  power  cannot 
be  delivered  through  the  conventional  ship's  service  power  distribution  system.  In  this  chapter, 
the  characteristics  of  a  shipboard  prime  power  distribution  system  capable  of  delivering  power 
to  these  FEL  components  are  explored. 

The  prime  power  supply  requirements  are  summarized  in  Table  8.  These  requirements 
are  established  to  be  consistent  with  the  FEL  power  requirements  of  Chapter  III. 


rf  sources 

dc  photocathode 

voltage 

100  kV 

500  kV 

power 

10  MW 

250  kW 

type 

dc 

dc 

duration 

15  sec 

15  sec 

Table  8  -  Prime  power  requirements  to  support  a  high-power  FEL 

The  value  of  100  kV  for  the  rf  source  input  voltage  was  chosen  to  support  the  high-power 
klystron  operation.  This  is  a  limiting  value  and  would  likely  be  about  40  kV  when  alternate  rf 
sources  are  developed. 

A.  POWER  GENERATION  AND  DISTRIBUTION 

The  FEL  system  will  require  a  dedicated  distribution  bus.  Therefore,  this  bus  is  not 
limited  to  the  normal  ship's  service  distribution  system  operating  parameters  of  450  V,  3 
phase,  and  60  Hz.  Several  studies  have  been  performed  to  outline  the  prime  power  system  for 
high-power,  electrically-powered  weapons  with  power  requirements  similar  to  FELs 
[52][53][54][55][56].  The  power  generation  and  distribution  aspects  of  these  studies  are  equally 
pertinent  to  this   discussion   and   are   summarized   in  this   section  for  continuity  purposes. 
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Follow-on  sections  will  address  components  specific  to  the  FEL  systems. 

The  prime  power  distribution  system  is  illustrated  in  Figure  16.  Power  generation  and 
distribution  component  weights  and  volume  are  summarized  in  Table  9.  Power  is  diverted  from 
a  propulsion  gas  turbine  engine  to  an  auxiliary  generator  to  power  the  FEL  distribution  bus. 
Modifications  to  the  pinion  that  connects  the  turbine  to  the  main  reduction  gears  have  been 
designed  so  that  the  turbine  power  is  automatically  diverted  to  the  auxiliary  generator  without 
operator  action  when  load  is  drawn  by  the  FEL  loads  [52] [53].  Simulation  studies  have  shown 
that  this  diversion  of  power  will  only  cause  a  momentary  reduction  of  ship  speed  of  a  few 
percent  when  operating  at  flank  speed  over  the  duration  of  the  laser  shot  [53].  The  IED 
generator  pictured  in  Figure  16  represents  the  electrical  generator  which  would  power  an 
electric  drive  ship.  This  element  would  be  deleted  for  a  ship  powered  with  mechanical  drive 
where  the  propeller  shaft  is  directly  coupled  to  the  reduction  gears. 

The  FEL  distribution  bus  would  operate  at  about  5  kV,  3  phase,  800  Hz.  As  bus 
frequency  is  increased,  the  size  and  weight  of  generator  and  transformer  components  are 
reduced.  However,  their  efficiency  is  also  reduced.  A  frequency  of  800  Hz  was  determined  to 
be  the  best  compromise  between  these  two  competing  factors  [52]  [55].  The  5  kV  bus  voltage 
was  chosen  due  to  conductor  sizing  limitations.  A  low  bus  voltage  causes  the  transmitted 
current  to  be  high  resulting  in  the  need  for  large  cable  conductors.  Similarly,  a  large  voltage 
results  in  increased  cable  size  due  to  insulation  requirements. 

B.  20  MVA  TRANSFORMER 

Computer  algorithms  have  been  developed  for  the  prediction  of  transformer  weights  for  a 
pulsed,  high-power  electrical  load.  Over  a  power  range  of  10-50  MW,  the  weight  of  an 
adiabatic  (no  forced  cooling)  transformer  is  given  by  the  expression  [57] 

Weight  =  0.0505  (r/120)  *  [0.693  +  0.307  (P/25)-019] 

*  [0.931  +  0.069  (Vo/100)13]  *  [0.242  +  0.758  /~°926]      (lbs/kW).     (4.1) 
For  a  1-5  MW  transformer,  the  weight  is  given  by 

Weight  =  0.1275  (//120)0281  *  [0.693  +  0.307  (P/25)-079] 

*  [0.931  +  0.069  (V^lOO)1 3]  *  [0.242  +  0.758  /"°-926]      (lbs/kW)      (4.2) 

where  r  is  the  total  run  time  (sec),  P  is  the  power  level  (MW),  V0  is  the  dc  load  voltage  (kV) 
assuming  rectification  without  phase-control  at  the  output  of  the  transformer,  and  /  is  the  bus 
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Figure  16   -  Prime  power  distribution  system.  After  Ref.  [52] 
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component 

weight 

volume 

(kg) 

(m3) 

gas  turbine 

20,400* 

46* 

Main  reduction  gears 

6800* 

3.8* 

auxiliary  generator 

8000 

4.3 

transmission  cable 

4500 

1.0 

auxiliary  support 

3300 

NA 

Table  9  -  Summary  of  power  generation  and  distribution  component  weights  and 
volume.  *  -  part  of  propulsion  system  and  does  not  add  to  FEL  system 
totals;   NA  -  not  available 


frequency  (kHz).  These  equations  assume  an  input  voltage  of  1-5  kV  and  have  been  validated 
for  output  voltages  of  20-200  kV,  run  times  of  1-120  seconds,  off  times  of  2-300  seconds,  and 
frequencies  up  to  2  kHz.  The  studies  which  produced  these  relations  also  showed  that  adiabatic 
transformers  provide  for  the  lowest  size  and  weight  transformer  characteristics  for  pulsed- 
power  applications. 

In  order  to  conform  with  Navy  standards,  all  components  must  be  derated.  The  derating 
factors  for  both  current  and  voltage  are  0.7.  The  transformer  which  will  supply  the  10  MW 
bus  for  the  rf  sources  was  derated  by  (0.7)2.  Therefore,  it  was  rated  at  20  MW.  From  equation 
4.1,  the  required  100  kV,  20  MW  transformer  would  weigh  620  lbs  (280  kg).  The  weight  of  a 
100  kV,  2  MW  transformer  would  be  157  lbs  (71  kg).  Therefore,  there  is  a  definite  weight 
advantage  to  using  a  single  transformer  to  supply  a  single  100  kV  bus  instead  of  using  a  few 
smaller  ones. 

A  volume  determination  can  be  made  from  the  method  described  in  Ref.  [55].  Figure  17 
in  Ref.  [55]  gives  the  volume  of  a  transformer  as  a  function  of  frequency  and  power  level. 
From  this  curve,  the  20  MW  transformer  would  occupy  approximately  5  m3. 

C.  100  KV  RECTIFIER 

The  rectifiers  were  assumed  to  be  composed  of  silicon-controlled  rectifiers  (SCRs)  instead 
of  diodes.    This  is  because  the  accelerator  rf  sources  and  dc  photocathode  voltages  need  to  be 
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controlled  to  tight  tolerances.  Use  of  a  phase-controlled  rectifier  allows  for  a  control  system  to 
regulate  the  output  voltage.  The  weight  of  the  100  kV  rectifier  was  estimated  to  be  1500  kg 
while  the  volume  was  estimated  to  be  8  m3.  The  analysis  used  to  determine  these  values  is 
presented  in  the  Appendix. 

D.  500  KV  POWER  SUPPLY 

The  500  kV  power  supply  must  provide  an  average  current  of  500  mA.  A  voltage 
multiplication  circuit  would  provide  better  performance  [12]  and  be  smaller  and  lighter  than  an 
equivalent  transformer  and  rectifier  combination.  The  weight  and  volume  of  a  500  kV  voltage 
multiplier  was  estimated  to  be  4100  kg  and  3  m3,  respectively.  The  analysis  used  to  determine 
these  values  is  presented  in  the  Appendix. 

E.  ENERGY  STORAGE 

An  important  consideration  for  any  weapons  system  is  how  it  operates  under  casualty 
conditions.  A  dedicated  energy  storage  element  could  provide  a  limited  capability  for  FEL 
operation  should  the  normal  power  distribution  be  interrupted.  Several  methods  of  energy 
storage  were  considered;  storage  batteries,  superconducting  magnetic  energy  storage,  flywheels, 
and  capacitor  banks.  Only  capacitor  banks  are  adequate  to  deliver  megawatts  of  power  within  a 
fraction  of  a  second  and  yet  be  compatible  with  the  shipboard  environment.  The  weight  and 
volume  of  high-energy  storage  capacitors  is  proportional  to  the  storage  capacity.  Modern 
capacitors  can  store  3  KJ/kg  [53]  and  it  is  anticipated  that  4  KJ/kg  is  achievable  with  some 
development  [54][58].  These  capacitors  would  occupy  a  volume  of  5  MJ/m3  [53][58]. 

The  run  time  of  15  seconds  for  the  FEL  is  based  on  the  time  available  to  engage  a  target 
and  not  on  the  laser  dwell  time  required  to  destroy  the  target.  Therefore,  it  may  not  be 
necessary  to  store  enough  energy  for  a  full  15  second  engagement.  20  MJ  of  stored  energy 
would  provide  for  a  limited  2  second  engagement.  This  amount  of  stored  energy  would 
require  a  capacitor  bank  that  weighs  5000  kg  and  occupies  4  m3.  This  capacitor  bank  would 
also  require  a  power  electronics  package  on  the  order  of  the  10  MW  rectifier  system  in  size 
and  weight.   The  detailed  estimates  of  this  package  were  not  made. 

F.  SUMMARY  OF  FEL  PRIME  POWER  DISTRIBUTION 

The  weight  and  volume  estimates  for  the  prime  power  distribution  system  is  summarized 
in  Table  10.  Those  components  which  are  normally  part  of  the  propulsion  train  are  not 
included.    The  transformer  and  rectifier  values  were  determined  for  a  100  kV  FEL  bus  voltage. 
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component 

weight                           volume 
(kg)                               (m3) 

power  generation 

transformer 

rectifier 

500  kV  power  supply 

15800                              5.3* 
280                                   5 
1500                                  8 

4100                                  3 

totals 

21680                              21.3 

Table  10  -  Prime  power  weight  and  volume  estimates.    *  -  does  not  include 
auxiliary  support  equipment 


These  weight  and  volume  estimates  can  be  reduced  appropriately  as  alternate  rf  sources 
become  available. 

Energy  storage  is  not  included  in  this  summary.  However,  its  inclusion  would  result  in 
an  increase  of  more  than  5000  kg  and  5  m3  to  the  totals.  This  would  be  a  significant  increase 
in  overall  system  size  and  weight. 
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V.  HIGH-POWER  FEL  OPERATION 

In  the  previous  chapters,  the  design  constraints  and  prime  power  considerations  for  a 
shipboard  high-power  FEL  have  been  outlined.  Various  electron  beam  parameters  have  been 
defined  so  that  an  FEL  could  provide  a  1  MW  optical  output.  In  this  chapter,  the  conversion 
of  electron  beam  energy  to  optical  energy  is  explored.  The  undulator  parameters  considered  in 
this  chapter  are  not  meant  to  specify  a  required  design  but  to  demonstrate  that  the  defined 
electron  beam  parameters  can,  in  fact,  produce  the  desired  output.  Both  the  undulator  and 
optical  resonator  parameters  are  also  constrained  by  the  overall  system  goals  of  minimum  size 
and  weight. 

The  electron  beam  parameters  that  must  be  considered  to  adequately  characterize  the  FEL 
operation  are  summarized  in  Table  1 1 . 


parameter 

value 

units 

Beam  energy,  Wb 

100 

MeV 

Peak  Current,  1^ 

360 

A 

M 

12 

- 

Micropulse  length,  Le  (400  MHz) 

42 

psec 

Micropulse  length,  Le  (200  MHz) 

84 

psec 

Beam  radius,  rb 

330 

urn 

Energy  spread,  Ay/y 

0.002 

- 

e  (rms) 

0.1 

mm-mrad 

Table  1 1  -  Electron  beam  parameters 

The  electron  beam  radius,  rb,  the  electron  beam  energy  spread,  Ay/y,  and  the  electron  beam 
emittance,  e,  are  estimates  made  from  existing  accelerator  electron  beams.  A  high  quality 
electron  beam  from  a  1.5  GHz  accelerator  with  a  peak  current  of  200  A  has  the  values 
rb  =  164  urn,  Ay/y  =  0.1%,  and  an  £  =  0.1  mm-mrad  [2].  The  election  beam  radius  and 
emittance  scale  approximately  as  the  current  [12].  It  has  already  been  established  that  emittance 
is  reduced  in  lower  rf  frequency  accelerators.  Furthermore,  electron  energy  spread  Ay  is 
conserved    as    an   electron   beam    is    accelerated    [22].    Therefore,    values   of   rb  =  330um, 
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Ay/y  =  0.002,  and  e  =  0.1  mm-mrad  are  reasonable  assumptions. 

A.  FEL  OPERATION  -  3  ^m  WAVELENGTH 

The  operating  wavelength  of  an  FEL  is  determined  by  the  resonance  condition,  equation 
2.6,  where  K  =  eB0k0/2m?ic2.  For  operation  at  X  =  3  p.m,  a  range  of  FEL  undulator 
parameters  are  summarized  in  Table  12.  These  parameters  were  defined  in  Chapter  II,  sections 
A  and  D. 


K 

K 

S 

Kl® 

B0  (peak) 

(cm) 

(kG) 

1 

4.7 

.48 

3.3 

71 

3 

2.6 

.43 

1.9 

13 

5 

1.9 

.39 

1.5 

5.7 

10 

1.1 

.28 

.96 

1.7 

Table  1 2  -  Undulator  parameters  for  FEL  operation  at  X  =  3  \un 


The  undulator  field  strength  is  limited  by  technology.  The  maximum  peak  undulator  field 
strength  that  can  be  produced  is  given  by  the  expression  [59] 

g 


(5.4" 


Bn  =  33.3  e 


1.8  —  ) 


(kG) 


(5.1) 


where  g  is  the  gap  between  the  pole  pieces  of  the  undulator  magnets.  The  strongest  field 
strength  that  has  been  manufactured  in  an  undulator  was  B0  =  18  kG  with  a  gfX0  =0.11  [60]. 
The  minimum  gap  distance  is  determined  by  the  transverse  dimension  of  the  optical  mode. 
Under  these  limitations  the  minimum  undulator  wavelength  that  can  be  considered  for  this  high 
power  FEL  is  X0  ~  3  cm  with  a  g  ~  5.4  mm. 

The  electron  micropulse  passing  through  the  undulator  is  several  hundred  slippage 
distances  long.  Hence,  the  optical  pulse  will  be  approximately  the  same  length  as  the 
micropulse  and  the  duty  cycle  between  pulses  will  also  be  equal.  For  the  parameters  of  Table 
11,  the  duty  cycle,  d,  will  be  about  1/720.  The  optical  pulse  amplitude  and  phase 
characteristics  will  be  periodic  along  the  length  of  the  pulse  with  a  period  of  about  one 
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slippage  distance.   Therefore,   simulation  of  the  FEL  electron-optical  interaction  over  one 
slippage  distance  is  sufficient  to  model  the  entire  pulse. 

The  relationship  between  average  power  output  of  the  optical  resonator  and  the  complex 
field  strength  is  given  by 


Pavg  =  nw0  d 


4n 


1         ~> 

ftnc  ~ 


4neNLKL  (£) 


(l-e-llQ)    \a\2     (kW) 


(5.2) 


where  Q  is  the  quality  factor  of  the  optical  resonator.  Over  the  length  of  the  undulator, 
however,  the  complex  field  strength  and  optical  mode  cross-sectional  area  are  not  constant  due 
to  diffraction  effects.  This  variation  is  a  direct  function  of  the  Rayleigh  length,  z0.  For  a 
high-power  FEL,  we  want  z0  and  w0  to  be  as  small  as  possible  to  minimize  the  power  density 
on  the  resonator  optics.  However,  if  w0  <  rb ,  then  the  gain  of  the  FEL  is  reduced  since  not  all 
of  the  electrons  can  interact  with  the  optical  mode.  As  explained  in  Chapter  III,  section  E,  this 
is  accounted  for  by  setting  j  =  j  F ,  where  F  =  (rhlw0)2,  when  simulating  the  wave  equation. 
The  value  vv0,  however,  is  dependent  on  the  Rayleigh  length  of  the  optical  resonator. 
Assuming  Gaussian  beam  propagation,  the  value  z0  »  L/VT2  minimizes  the  volume  for  the 
optical  mode  over  the  length  of  the  undulator  and  provides  for  maximum  gain  when  w(J  >  rh 
[5].   With  z0  =  L/VT2,  we  can  use  equation  5.2  and  neglect  the  variation  in  la  I. 

The  optical  parameters  for  several  undulators  are  listed  in  Table  13.   The  Rayleigh  length 
of  the  resonator  is  assumed  to  be  z0  =  L/VT2. 
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3 

15 

.45 

.35 

6.3 

25 

.75 

.46 

17.6 

4 

15 

.6 

.41 

6.3 

25 

1 

.53 

17.4 

5 

15 

.75 

.46 

6.2 

25 

1.25 

.59 

17.1 

Table  13  -  Optical  parameters  for  FEL  operation  at  X  =  3  [im 
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Several  simulations  were  performed  to  determine  the  relationship  between  j,  la  I,  and  Q .  The 
results  of  these  simulations  for  values  of  ;  =  4,  7,  and  13  are  shown  as  curves  in  Figure  17. 
Overlayed  with  these  curves  is  the  \a  I  vs.  Q  relation  determined  by  equation  5.2  for  a  1  MW 
optical  output  and  an  undulator  with  X0  =  3  cm,  N  =  16,  and  j  =  7. 
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Figure  17  -  la  I  vs.  Q  dependence  for  several  values  of  j.  The  curve  marked 
"P  =  1  MW"  is  the  Ifl  I  vs.  Q  relationship  for  an  undulator  with  j  =  7,  K0  =  3cm , 
and  N  =  16.  The  intersection  of  the  two  solid  curves  is  the  required  operating 
point 


The  optical  mode  within  the  resonator  must  have  an  la  I  «  17  and  the  resonator  must  have  an 
outcoupling  of  20%  (Q  ~  4.5).  From  Figure  17,  it  is  obvious  that  more  than  1  MW  can  be 
extracted  from  the  electron  beam  by  increasing  the  resonator  Q.  However,  this  would  result  in 
an  undulator  extraction  efficiency  of  more  than  2%. 

The  average  power  within  the  resonator  cavity  is  5  MW  and  is  contained  within  an 
optical  mode  size  of  wn  =  0.364  mm.  Assuming  a  Gaussian  profile,  essentially  all  of  the 
power  (<    100  W)  is  contained  within  a  cross-section  of  8u0  =  2.7  mm  at  the  ends  of  the 
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undulator  where  the  spot  size  is  largest.  Allowing  for  alignment  tolerances,  the  undulator  gap 
distance  g  =  5.4  mm  is  adequate. 

B.  ANALYSIS  OF  FEL  OPERATION 

The  FEL  saturates  at  an  average  power  of  1  MW  for  j  =  1  and  Q  ~  4.5.  The  steady 
state  characteristics  of  the  FEL  interaction  are  shown  in  Figures  17,  18,  and  19.  Figure  18  is 
the  gain  and  phase  curve,  as  defined  in  section  H  of  Chapter  II,  for  this  FEL  near  saturation. 
The  undulator's  sensitivity  to  electron  beam  quality,  as  explained  in  section  G  of  Chapter  II,  is 
proportional  to  the  number  of  periods.  Therefore,  the  gain  curve  is  virtually  insensitive  to 
beam  quality  since  the  undulator  only  has  16  periods.  Figure  19  is  the  phase  space  evolution  of 
an  electron  pulse  at  saturation.  Approximately  60%  of  the  electrons  are  trapped  in  closed 
orbits.  The  trapped  electrons  leave  the  undulator  with  phase  velocities  concentrated  about  v  ~ 
-5  while  the  non-trapped  electrons  exit  with  phase  velocities  spread  about  v  ~  5. 
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Figure  18   -   Gain  spectrum  at  saturation  for  operation  at  A,  =  3  (im 
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Figure  19  -   Phase  space  evolution  at  saturation  for  operation  at  A,  =  3  (im 


The  startup  characteristics  of  the  FEL  are  illustrated  in  Figure  20.  The  evolution  of  the 
optical  mode  over  n  =  150  micropulses  beginning  from  spontaneous  emission  is  shown.  The 
left  window  shows  the  complex  field  strength  over  one  slippage  distance,  the  center  window 
shows  the  power  spectrum,  and  the  right  window  shows  the  electron  distribution.  The  electron 
distribution  of  this  simulation  agrees  with  the  phase  space  distribution  of  Figure  19  in  that  the 
electrons  are  concentrated  in  two  distinct  phase  velocity  groups.  The  power  spectrum  has  a 
single  mode  and  so  there  is  no  modulation  of  the  complex  field  strength  envelope.    The 


steady-state  complex  field  strength  is  la  I  =  xP(n)  ~  17.  The  optical  field  reaches  saturation  in 
about  50  passes  (1.5  p,sec  for  an  rf  frequency  of  400  MHz  or  3  |isec  for  an  if  frequency  of 
200  MHz)  for  the  FEL  to  reach  the  steady-state  power  level  of  1  MW.  In  a  combat 
environment,  this  is  essentially  an  instantaneous  engagement. 
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Figure  20  -  Optical  field  characteristics  at  saturation  for  operation  at  A,  =  3  (im 

C.  FEL  OPERATION  -  10  (im  WAVELENGTH 

In  the  maritime  environment,  there  are  two  principal  wavelength  bands  in  the  infrared 
(IR)  suitable  for  laser  propagation.  These  are  the  3-5  p.m  short  wave  infrared  (S WTR)  and  the 
8-12  p.m  long  wave  infrared  (LWTR)  bands.  Now  that  operation  at  3  \un  with  a  1  MW  output 
has  been  demonstrated,  operation  in  the  LWTR  band  at  a  wavelength  about  10  p.m,  using  the 
same  undulator,  is  explored. 

To  maintain  resonance  at  A.  =  10  (im,  the  beam  energy  must  be  reduced  to  54.8  MeV. 
The  dimensionless  current  and  optical  mode  waist  are  now  y  =  13  and  w  =  .66  mm, 
respectively.  If  we  lay  a  curve  corresponding  to  a  1  MW  output  with  the  new  parameters  on 
top  of  Figure  17,  we  find  that  \a  I  ~  30  and  Q  ~  4.5.    Therefore,  if  the  resonator  optics  can  be 
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made  to  have  the  same  outcoupling  at  the  LWIR  wavelengths,  the  same  undulator  and 
resonator  cavity  could  be  used  to  operate  in  both  bands. 

The  simulation  results  equivalent  to  Figures  18,  19,  and  20  reveal  only  a  few 
characteristics  different  than  for  the  X  =  3|i//z  case.  There  are  four  principal  longitudinal 
modes  excited  by  the  micropulses.  Therefore,  there  is  a  modulation  on  the  complex  field 
envelope.  This  results  in  a  startup  time  increase  for  the  FEL  of  a  factor  of  three.  These 
features,  while  interesting,  do  not  affect  the  engineering  and  feasibility  results  of  concern  in 
this  analysis. 

Despite  the  apparent  wide  bandwidth  in  operating  wavelengths,  there  are  many 
repercussions  which  must  be  considered.  First,  to  reduce  the  beam  energy  by  about  one-half, 
the  accelerator  E  is  reduced  in  half.  Therefore,  the  energy  storage  in  each  accelerator  cavity  is 
reduced  by  a  factor  of  four.  This  means  that  AUIU  increases  by  a  factor  of  four  resulting  in 
an  equal  fractional  increase  in  the  rf  field  distortions  caused  by  an  accelerating  micropulse. 
This  will  likely  be  unsatisfactory  for  proper  beam  transport  in  RT  cavities.  SC  cavities  will  be 
more  resistant  to  this  effect. 

A  second  consideration  is  the  optical  mode  size  with  respect  to  the  undulator  gap 
distance.  In  the  event  of  perfect  alignment,  the  undulator  will  not  penetrate  into  the  Gaussian 
profile  of  the  optical  energy.  However,  if  the  alignment  is  off  by  as  little  as  1  mm,  the 
undulator  end  magnets  will  penetrate  into  space  normally  occupied  by  several  kilowatts  of 
power.  This  can  be  compensated  for  by  increasing  the  undulator  gap  distance.  However,  as  a 
result  of  the  corresponding  changes  in  the  undulator  parameters,  the  undulator  will  be  longer, 
have  an  increased  sensitivity  to  electron  beam  quality,  and  the  Q  value  of  the  optical  resonator 
will  increase.  This  will  increase  the  power  density  incident  on  the  mirrors  for  the  same  mirror 
separation.  Therefore,  if  operation  in  the  LWIR  is  desired,  this  tolerance  must  be  considered 
when  finalizing  the  undulator  and  resonator  design. 

D.  RESONATOR  OPTICS 

One  of  the  principal  disadvantages  of  the  FEL  for  use  as  a  laser  weapon  is  that  the 
optical  waist  is  very  small.  Therefore,  the  power  density  incident  on  the  resonator  optics  is 
large.  Figure  21  shows  the  average  power  density  incident  on  the  resonator  mirrors  as  a 
function  of  mirror  separation  for  the  undulator  with  X0  =  3  cm  and  N  =  16.  With  an  18  meter 
mirror  separation,  the  mirror  power  density  is  70  kw/cm 2  at  X  =  3  u.m  and  20  kW/cm :  at 
A.  =  10|im.  Modern  optical  components  can  support  a  power  density  in  excess  of  100 
kW/cm2  [51]. 
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Figure  21  -  Power  density  on  the  resonator  mirrors  for  operation  at  A.  =  3  urn  and 
X  =  10  |im  (  *  -  denotes  aLlowable  mirror  placements  consistent  with  the  rf 
frequency  choices  in  Table  7) 


In  the  shipboard  environment,  space  and  volume  allocation  is  of  particular  concern  to  the 
ship  design  engineers.  The  mirror  power  density  limitations  require  that  the  mirrors  be  about 
18  m  apart.  However,  it  may  be  possible  to  reduce  the  resonator  length  with  mirrors  that 
exceed  the  100  kW/cm2  performance  level.  While  this  is  a  considerable  length  within  the  ship, 
the  actual  volume  of  the  resonator  is  quite  small  since  the  optical  mode  can  be  contained 
within  a  cross-section  of  about  20  cm.  The  hardware  for  the  active  stabilization  of  the  optics  to 
maintain  resonator  alignment  during  ship  vibrations,  flexing,  and  torsional  rotations  will 
constitute  the  bulk  of  the  resonator  volume.  Therefore,  the  only  significant  amount  of  space 
that  will  be  needed  will  be  at  the  mirror  locations  themselves. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  FEL  is  a  promising  candidate  for  a  possible  high-average  power  shipboard  laser 
weapon.  Unlike  other  candidates,  such  as  chemical  lasers,  the  FEL  is  electrically  powered. 
An  additional  electrical  power  distribution  bus  would  be  required  aboard  a  ship  if  an  FEL 
weapon  was  to  be  deployed. 

An  FEL  with  a  1  MW  laser  output  would  require  6-10  MW  of  electrical  power  from  the 
prime  power  distribution  system.  A  gas  turbine  naval  engineering  plant  could  be  modified  to 
supply  this  power.  Power  from  a  propulsion  turbine  could  be  diverted  to  power  the  dedicated 
FEL  distribution  bus.  The  added  distribution  system  would  increase  the  ship's  displacement  by 
22  tons  and  require  22  m3  of  the  ship's  volume  to  implement  The  power  supply  that  powers 
the  if  sources  for  the  linear  accelerator  were  rated  at  an  output  voltage  of  100  kV  and  an 
output  power  of  10  MW.  This  power  supply  is  only  82%  efficient  and  requires  a  large 
external  cooling  capability.  The  efficiency  can  be  improved,  however,  with  the  development  of 
alternative  rf  sources  which  operate  at  half  of  this  voltage.  The  development  of  these  alternate 
if  power  supplies  is  recommended. 

The  linear  accelerator  driving  the  FEL  would  employ  energy  recovery  to  increase  the 
power  efficiency  and  reduce  the  operating  radiation  levels.  Consequently,  the  electron  beam 
has  a  large  average  current.  Furthermore,  the  accelerator  would  operate  at  lower  if  frequencies 
than  most  linear  accelerators.  Therefore,  the  electron  pulse  lengths  are  longer  than  those 
commonly  used  in  most  FELs.  The  high-average  current  and  long  pulse  lengths  reduce  the 
number  of  undulator  periods  required  to  extract  the  1  MW  output  power  from  the  electron 
beam.   Therefore,  the  FEL  has  a  low  sensitivity  to  electron  beam  quality. 

The  FEL  examined  in  this  thesis  requires  further  development  before  it  can  be  built. 
Some  of  the  technologies  that  require  further  development  are:  energy  recovery  employment, 
control  system  requirements,  and  the  accelerator  electron  injector  output  capacity.  Additional 
experiments  are  needed  to  determine  the  limitations  of  energy  recovery  including  if  phase  and 
amplitude  stability  requirements  as  a  function  of  FEL  undulator  energy  extraction  efficiency. 
Control  systems  are  required  to  regulate  the  input  power  to  if  sources  and  the  rf  input  to 
accelerator  components.  Additional  work  is  necessary  to  adequately  control  the  electron  beam 
transport  instabilities  that  will  result  from  high-average  currents.  Finally,  electron  injector 
technology  can  currently  produce  average  currents  of  tens  of  milliamps.  Advancements  are 
needed  for  the  shipboard  FEL  weapon.  Further  development  of  these  FEL  technologies  is 
recommended. 
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APPENDIX  -  POWER  SUPPLY  CALCULATIONS 

This  appendix  outlines  the  procedure  used  to  estimate  the  weight  and  volume  of  the  FEL 
power  supplies  in  Chapter  IV.  Only  commercially  available,  off-the-shelf  components  were 
considered  in  the  analysis.  No  attempt  was  made  to  determine  the  savings  that  may  be 
possible  with  a  product  development  program.  Therefore,  the  results  are  a  conservative 
estimate.  The  components  selected  in  the  analysis  are  meant  to  be  representative  components. 
It  should  not  be  assumed  that  these  components  are  recommended  for  procurement.  The 
voltage  and  current  ratings  for  all  of  the  components  used  in  the  analysis  were  derated  by  the 
Navy  derating  factor  of  0.7.   The  dc  link  and  filter  components  were  not  considered. 

A.  100  KV,  3  PHASE,  PHASE-CONTROLLED  RECTIFIER 

The  functional  circuit  diagram  for  a  3  phase,  phase-controlled  rectifier  is  shown  in  Figure 
22. 
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Figure  22   -   Functional  circuit  diagram  for  a  3  phase,  phase-controlled  rectifier 

Commercially  available  silicon  controlled  rectifiers  (SCRs)  are  not  rated  to  support  a  100  kV 
blocking  voltage.  Several  SCRs  must  be  connected  in  series  to  accomplish  the  function  of 
each  SCR  in  Figure  22.  Each  of  these  series-connected  SCRs  must  be  protected  with  a 
snubber  circuit  and  balancing  resistor.  The  series  connection  of  the  SCRs  is  illustrated  in 
Figure  23.  Forced -com  mutation  circuits  for  the  SCRs  are  not  included.    It  is  assumed  that  the 
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Figure  23   -   Series-connected  SCRs  with  protective  balancing  resistors  and  snubbers 

blocking  voltage  is  sufficient  to  rapidly  reverse  bias  the  SCRs. 

1.  SCR  Selection 

The  rectifier  has  an  output  dc  voltage,  Vdc ,  of 


Vdc  = 


Vm  cosa=  100  kV 


(A.l) 


where  Vm  is  the  peak  line-to-line  input  ac  voltage  and  a  is  the  SCR  firing  angle.  Since  the 
output  dc  voltage  is  reduced  by  only  13.4%  when  a  is  adjusted  from  0-30°,  Vm  =  105  kV  is 
chosen  so  that  Vdc  =  100  kV  for  a  -  0°.  To  supply  10  MW  of  power,  the  output  dc  current, 
Idc,  must  be  100  A.  The  peak  reverse  blocking  voltage  for  the  SCR,  Vrr,  is  Vm  and  the 
conduction  current,  /,  is  100  A.  The  derated  values  of  Vrr  and  /  are  Vdrr  =  150  kV  and  ld  = 
143  A,  respectively. 

Using  Ref.  [61],  the  SCR  model  T8K7  was  chosen  for  the  SCR  component.  The  ratings 
of  this  SCR  are  listed  in  Table  14.  The  value  VRRM  is  the  repetitive  reverse  blocking  voltage, 
h(av}  ^  the  P&k  average  conduction  current,  lGT  is  the  gate  current,  IRRM  is  the  leakage 
current,  /^^  is  the  peak  repetitive  surge  current,  and  id  is  the  turn-on  delay  time  of  the  SCR. 
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Vrrm 
(kV) 

'T(av) 

(A) 

(mA) 

Irrm 
(mA) 

4.5 

350 

200 

75 

hsM 
(A) 

u.sec 

Weight 
(kg) 

Volume 
(cm3) 

1500 

2 

0.227 

72.9 

Table  14  -   Rating  values  for  SCR  model  T8K7 

Allowing  for  the  failure  of  one  SCR,  35  SCRs  are  required  to  construct  each  of  the  six 
branches  in  Figure  22.  Therefore,  a  total  of  210  SCRs  are  needed.  The  total  SCR  conduction 
losses  are  19  kW. 


2.  Balancing  Resistors 

Balancing  resistors  are  required  to  ensure  that  all  series-connected  SCRs  block  the  same 
voltage.  Otherwise,  variations  in  leakage  currents  could  cause  an  excessive  voltage  to  build  up 
across  one  SCR.  The  actual  blocking  voltage  for  each  SCR  is  Vblock  =  3.1  kV.  The  value  of 
Rh  is  calculated  as 


Rh  = 


block 


The  required  power  rating  of  this  resistor  is 


=  4  kQ 


(A.2) 


/V  = 


-  =  4.8kW 

*"         Rh 


(A.3) 


where  Vd  =  4.5  kV  is  the  derated  voltage  rating  of  the  SCRs.  Multiple  resistors  will  have  to 
be  placed  in  series  to  meet  this  power  rating.  From  Ref  [62],  a  suitable  resistor  is  the  type 
270.    The  ratings  of  this  resistor  are  listed  in  Table  15.  The  actual  power  dissipation  of  this 
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Table  15  -  Rating  values  for  the  type  270  resistor 

resistor  is  VtiocklRb  =  2.4kW .  One  balancing  resistor  is  needed  for  each  series-connected  SCR 
and  five  type  270  resistors  are  required  per  balancing  resistor.  Therefore,  a  total  of  1050 
resistors  are  needed  with  a  total  power  dissipation  of  500  kW. 

3.  Snubber  Components 

In  this  rectifier,  the  snubber  provides  voltage  and  current  protection  for  series-connected 
SCRs  during  firing  and  turn-off.  If  one  or  more  SCRs  latch  slower  than  the  rest  of  the  SCRs 
in  a  branch,  an  excessive  voltage  could  be  applied  across  these  SCRs.  The  snubber  capacitor 
provides  the  charge  to  the  conducting  SCRs  during  the  period  when  these  slower  SCRs  still 
have  not  latched.  When  the  SCRs  are  first  conducting,  though,  the  capacitor  would  be 
effectively  shorted  causing  a  large  current  surge.  The  snubber  resistor  is  placed  in  series  with 
the  capacitor  to  limit  the  discharge. 

There  is  no  definitive  way  to  analytically  determine  rectifier  snubber  component  values. 
However,  a  conservative  methodology  can  be  used  to  estimate  these  values  with  a  design 
margin.   The  following  method  describes  how  the  snubber  component  values  were  determined. 

The  latching  delay  time  for  the  T8K7  SCR  is  2  (isec.  Assuming  the  snubber  capacitor 
must  supply  the  full  conduction,  the  minimum  capacitance  necessary  to  store  the  required 
charge  is 

Csjnin   ="^=0.09  UP  (A.4) 

The  snubber  resistor  must  limit  the  capacitor  current  to  within  the  surge  current  rating  of  the 
SCR.  The  derated  surge  current  rating  of  the  SCR,  Idsm ,  is  1050  A.  The  minimum  snubber 
resistance  is  calculated  by 
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R. 


Idsm    ~  I d 


=  5Q 


(A.5) 


Assuming  that  the  SCR  switches  instantaneously,  the  power  rating  of  the  snubber  resistor  is 
required  to  be 


Pr.   ~  Irms  Rs    ~ 


V} 


R. 


RsCsfa 


(A.6) 


where  fac  is  the  rectifier  input  ac  frequency.  Equations  A.4,  A.5,  and  A.6  include  a  margin  of 
conservatism.  The  resistor  power  losses  are  independent  of  the  resistance  value  but  are  directly 
proportional  to  the  ac  bus  frequency.  To  minimize  the  power  losses  while  maintaining  the 
effectiveness  of  the  snubber,  the  values  of  Rs  =  5  Q  and  Cs  =  0.2  ^F  are  required.  The 
required  voltage  rating  of  the  capacitor  is  4.5  kV,  and  the  resistor  is  to  have  a  power  rating  of 
3.2  kW.  The  value  fac  =  800  Hz  is  assumed. 

From  Ref.  [63],  a  suitable  capacitor  design  was  selected.   The  capacitor  ratings  are  listed 
in  Table  16. 


Voltage 
(kV) 


Capacitance 
(UF) 


Weight 
(kg) 


Volume 

(cm3) 


Table  16  -   Rating  values  for  a  0.2  \iF,  4.5  kV  capacitor 

There  are  210  required  snubber  capacitors.   Therefore,  a  total  of  210  capacitors  are  required. 

The  type  270  was  selected  for  the  snubber  resistor.  Four  type  270  resistors  are  needed 
for  each  snubber.  Therefore,  a  total  of  840  resistors  are  required.  The  power  loss  of  each 
snubber  resistor  is  V^lockCsfac  =  1.54kW,  so  that  the  total  power  dissipation  by  all  of  the 
snubbers  is  1290  kW. 
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4.  SCR  Gating  Power  Supply 

A  current  pulse  applied  to  the  gate  of  an  SCR  is  required  in  order  for  it  to  fire.  The 
T8K7  SCR  requires  a  current  pulse  of  200  mA  and  dissipates  5  W  at  the  gate  in  order  to  fire. 
Since  there  are  35  SCRs  in  each  of  the  six  rectifier  branches,  a  total  of  7  A  and  175  W  is 
required  to  fire  all  of  the  SCRs  in  a  branch.  The  derated  values  of  required  gate  power  supply 
current  and  power  are  10  A  and  250  W,  respectively. 

From  Ref.  [64],  Acopian  power  supply  model  A28H1400  was  selected.  The  ratings  of 
this  power  supply  are  listed  in  Table  17. 


Voltage 
(kV) 

Current 
(A) 

Power 
(W) 

Weight 

(kg) 

Volume 
(cm3) 

28 

12 

336 

11.8 

9900 

Table  17  -   Rating  values  for  Acopian  power  supply  model  A28H1400 


A  total  of  six  power  supplies  would  be  required. 


B.  VOLTAGE  MULTIPLICATION  CIRCUIT 

The  fundamental  circuit  diagram  for  a  500  kV  voltage  multiplier  is  shown  in  Figure  24. 
Figure  24  is  an  8x  voltage  multiplier.  An  increased  multiplication  factor,  Fv,  is  obtained  by 
cascading  additional  stages.  To  obtain  a  500  kV  output  with  an  input  line  voltage  of  5  kV 
(rms),  a  multiplication  factor  of  Fv  =  72  is  required.  It  is  important  to  know  that  if  any  of 
these  stages  fail,  the  voltage  multiplication  factor  is  reduced  by  two.  The  next  stage  in  the 
cascade  must  be  rated  to  assume  the  load  of  the  failed  stage  or  it  may  also  fail.  Therefore, 
each  diode  and  capacitor  must  be  designed  to  the  limiting  position  in  the  circuit.  The  diodes 
must  be  designed  for  a  reverse  blocking  voltage  of  2Ep  and  a  current  of  FvIdc .  The  capacitors 
must  be  rated  for  a  voltage  of  3Ep .  The  derated  peak  surge  current,  to  which  all  of  the  circuit 
components  must  be  rated,  is  given  by 


1dm  =^facCEp  =228.5/1. 


(A.9) 
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Figure  24   -   Functional  diagram  for  a  voltage  multiplier  circuit 

1.  Capacitor  selection 

The  effective  output  capacitance,  Cfff ,  of  the  voltage  multiplier  needed  to  supply  a 
current  Idc  is  given  by 


c  f* 

where  AV0  is  the  allowed  output  voltage  ripple.   The  value  Ccff  of  the  circuit  is 

r//  FJFv-2) 


(A.7) 


(A.8) 


where  C  is  the  capacitance  of  each  capacitor.  Given  values  of  fac  =  800  Hz,  ldc  =  500  mA, 
Fv  =  72,  and  arbitrarily  choosing  AV0  =  5  kV,  then  each  capacitor  in  Figure  24  must  have  a 
capacitance  C  =  2.25  |iF.  The  derated  voltage  rating  of  the  capacitors  is  30.3  kV.  The  total 
number  of  these  capacitors  is  1.5FV  -  2  =  106  capacitors. 

Using  Ref.  [63],  capacitor  model  30F1377  was  selected  for  the  capacitor  component. 
The  ratings  of  this  capacitor  are  listed  in  Table  18. 
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voltage 
(kV) 


Capacitance 
(J1F) 


Weight 
(kg) 


Volume 
(cm3) 


30 


2.37 


19 


11500 


Table  18  -   Rating  values  for  capacitor  model  30F1377 

A  total  of  106  capacitors  is  required. 

2.  Diode  selection 

Each  diode  in  Figure  24  must  be  rated  for  the  limiting  current  and  voltage  values  in  the 
circuit.  The  required  diode  ratings  are  Id  =  FvIdc/0.1  =  51.4  A,  Vdrr  =  20.2  kV,  and  Idsm  = 
228.5  A.  Using  Ref.  [61],  diode  R4051070  was  selected.  The  ratings  of  this  diode  are  listed 
in  Table  19. 


Table  19  -   Rating  values  for  diode  R4051070 

It  is  necessary  to  place  several  of  these  diodes  in  series  to  accomplish  the  function  of  each 
diode  in  Figure  24.  Allowing  for  a  single  failure  within  the  diode  branch,  22  diodes  are 
required  per  branch.  There  are  Fv  diode  branches  so  that  a  total  of  1584  R4051070  diodes  are 
required.   The  total  conduction  losses  are  50  kW. 
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3.  Balancing  resistor  selection 

The  peak  blocking  voltage  of  each  diode  is  Vblock  =  2Ep/22  =  643  V.  Using  equation 
A.2,  the  required  value  for  the  balancing  resistor  is  Rb  =  643  kQ.  The  required  power  rating 
for  this  resistor  is  V}lRb  =  1.3  W.  Using  Ref.  [65],  resistor  620KW-2-ND  was  selected.  The 
ratings  of  this  resistor  are  listed  in  Table  20. 


Table  20  -   Rating  values  for  resistor  620KW-2-ND 

One  of  these  resistors  is  required  for  each  of  the  R4051070  diodes.  Each  of  these  resistors 
dissipates  Vb2hck/Rh  =0.64W.  Therefore,  1584  resistors  are  required  with  a  total  power 
dissipation  of  1  kW. 

C.  SUMMARY 

The  total  weight,  volume,  and  power  losses  for  the  100  kV  rectifier  and  500  kV  voltage 
multiplier  power  supplies  are  listed  in  Table  21. 


100  kV 

500  kV 

units 

Weight 

1500 

4100 

kg 

Volume 

8 

3 

m> 

power  losses 

1800 

50 

kW 

Table  21  -   Power  supply  summary 
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The  weight  estimate  was  made  by  multiplying  the  component  totals  by  a  factor  of  two. 
Power  supplies  that  operate  above  about  40  kV  are  required  to  have  enclosures  which  are  filled 
with  high  voltage  oil,  SF6,  or  some  other  means  of  protection  between  components.  Due  to  the 
large  amount  of  heat  losses,  the  10  MW  power  supply  is  assumed  to  be  filled  with  oil  which 
can  also  be  used  as  a  cooling  medium.  The  500  kV  power  supply  will  need  high  voltage 
shielding  to  isolate  it  from  the  environment.  The  volume  estimates  include  a  stacking  factor  of 
2.5.  This  factor  accounts  for  the  high  voltage  shielding  and  unused  space  within  the  power 
supply  after  the  components  are  assembled. 

The  power  efficiency  of  the  10  MW  power  supply  is  about  82%.  While  this  is 
acceptable  for  desktop  size  electronics,  the  1.8  MW  of  power  losses  place  an  added  design 
burden  of  forced  cooling  for  the  power  supply.  The  efficiency  of  the  power  supply  can  be 
increased  significantly  with  the  development  of  alternative  accelerator  rf  sources,  as  described 
in  Chapter  III,  which  operate  at  substantially  lower  voltages. 
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